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LIGHTLY  ARMORED  STRUCTURE  VULNERABILITY 
ESTIMATION  METHODOLOGY  (LASVEM) 


\ 


1.  INTRODUCTION 


The  USAMSAA  has  developed  a computer  model  which  assists  in  the 
evaluation  of  lightly  armored  structures,  such  as  armored  personnel 
carriers,  by  computing  various  measures  of  vulnerability  of  those  structures 
when  subjected  to  indirect  fire,  high  explosive  (not  ICM)  field  artillery 
projectiles.  The  basic  concept  is  that  the  primary  function  of  such  a 
structure,  with  respect  to  field  artillery  fire, -is  to  prevent  the  intrusion 
of  fragments  into  an  enclosure;  therefore,  a fundamental  measure  of  such 
a structure's  vulnerability  to  artillery  fire  is  the  number  of  perforations 
which  it  sustains  when  subjected  to  such  fires.  Such  a quantity,  however, 
may  not  be  especially  useful  to  decision  makers  at  various  levels.  Hence, 
a family- of  statistical  parameters,  all  based  upon  the  expected  number  of 
perforations,  has  been  developed  to  address  the  vulnerability  of  a group  of 
structures  to  various  volumes  of  tactically  delivered  field  artillery  fires. 


The  program  was  originally  developed  for  use  on  the  USA  Ballistic 
Research  Laboratory  (USABRL)  BRLESC  I/II  computers.  The  characteristics  of 
those  computers  resulted  in  restrictions  on  the  size  of  certain  matrices 
used  in  the  program  and  in  the  use  of  a magnetic  tape  oriented  file  system. 
As  documented  herein,  the  program  has  been  modified  to  run  on  the  new 
USABRL  CYBER  173/76  computer.  The  file  system  is  now  random  access,  but 
the  size  of  the  above  mentioned  matrices  has  not  been  increased  in  the 
desire  to  minimize  the  changes  made  to  an  already  debugged  program.  This 
is  the  source  of  certain  remarks  in  this  report  concerning  restraints  due 
to  machine  size. 

This  document  is  intended  to  serve  as  a user's  guide  so  that 
other  organizations  may  achieve  a capability  in  this  area.  It  is  also 
intended  to  serve  as  a programmer's  guide  for  those  changes  to  the  pro- 
gram which  may  be  necessary  or  desired. 

2.  OVERVIEW  OF  THE  COMPUTATIONAL  TECHNIQUE 


2.1  The  Program. 

The  program  is  divided  into  four  main  sections.  Section  I uses 
the  methodology  of  Thor  reports  47  and  70  (references  1 and  2)  to  produce 
a description  of  the  damage  sustained  by  a single  plate  or  plate  pair  from 
the  detonation  of  a projectile  at  a fixed  point  in  space  with  a fixed 
angle  of  fall  and  terminal  velocity,  for  all  locations  of  the  plate  or 
plate  pair.  Section  II  combines  the  descriptions  of  damage  sustained  by 
several  plates  (output  of  Section  I)  to  obtain  a description  of  damage 
sustained  by  a complete  structure.  Section  III  uses  the  descriptions  of 
damage  sustained  by  a complete  structure  from  the  detonation  of  a projec- 
tile to  obtain  various  statistical  measures  of  the  damage  sustained  by 
specified  arrays  of  target  structures  when  these  arrays  are  subjected  to 
tactically  delivered  field  artillery  fire. 
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Whereas  the  first  three  sections  of  the  program  deal  with  the 
number  of  perforations  sustained,  Section  IV  assesses,  in  a rather 
simplistic  manner,  the  effect  of  a perforation  upon  people  within  the 
enclosure.  Additionally,  it  determines,  without  reference  to  the  detailed 
and  time  consuming  computations  of  Section  I,  whether  any  of  the  fragments 
have  the  capability  of  ever  penetrating  (except  by  direct  hit  or  very  near 
miss  of  the  projectile)  the  plate  under  consideration. 

2.2  Damage  Criteria. 

The  function  of  a lightly  armored  structure  with  respect  to  field 
artillery  fire  is  to  prevent  the  intrusion  of  fragments  into  an  enclosure. 
(Protection  against  direct  hits  and  very  near  misses  seems  impossible  to 
achieve  with  current  light  armors.)  Therefore,  the  damage  criteria  used 
herein  are  always  of  the  form,  "sustain  at  least  N perforations",  and  the 
probability  statements  obtained  are  always  of  the  form,  "probability  that  M 
targets  each  sustain  at  least  N perforations".  For  brevity,  this  is  often 
expressed  as,  "probability  that  M vehicles  are  damaged",  where  the  damage 
criterion  has  been  specified  as  above. 

A minor  problem  arises  with  damage  criteria  of  this  form  if 
one  wishes  to  differentiate  between  the  damage  (which  may  be  major)  caused 
by  individual  fragments  and  the  truly  catastrophic  damage  caused  by  direct 
hits  or  very  near  misses.  Within  this  model,  the  problem  has  been  elim- 
inated by  forcing  a direct  hit  or  very  near  miss  to  produce  a very  large 
expected  number  of  perforations  (currently  1(T).  Thus,  direct  hits  and  very 
near  misses  are  addressed  within  the  established  framework,  and  yet  are 
distinguished  merely  by  using  a criterion  "sustain  at  least  N perforations", 
where  N is  large  but  still  small  with  respect  to  10  . 

2.3  Measures  of  Vulnerability. 


Sections  I and  II  produce  an  expected  perforations  matrix  for  a 
single  target  subjected  to  the  fragmentation  of  one  projectile  at  a fixed 
point  in  space  with  a fixed  angle  of  fall  and  terminal  velocity.  Each 
element  of  this  matrix  is  the  expected  number  of  perforations  sustained  by 
the  target  when  it  is  located  at  a point  in  the  ground  plane  corresponding 
to  the  matrix  coordinates,  (i,j). 

The  probability  of  damage  matrix  is  obtained  from  the  expected 
perforations  matrix.  With  a damage  criterion  of  the  form,  "sustain  at 
least  N perforations",  and  the  assumption  that  a Poisson  distribution 
applies  for  each  location  of  the  target,  each  element  of  this  matrix  is  the 
probability  that  the  target  sustains  damage  at  that  location.  There  is  a 
one-to-one  correspondence  between  elements  of  the  expected  perforations 
matrix  and  elements  of  the  probability  of  damage  matrix. 


From  the  probability  of  damage  matrix  are  obtained: 

a.  The  lethal  area  (reference  3). 

b.  The  average  probability  of  damage  within  the  area  over  which 
the  probability  of  damage  is  non-zero. 

c.  The  area  over  which  the  probability  of  damage  is  non-zero. 

d.  A damage  function. 

Section  III  uses  the  methodology  of  reference  4 to  address  the 
vulnerability  of  an  array  of  targets  to  various  volumes  of  artillery  fire. 
Using  damage  functions  (one  for  each  HOB  in  the  HOB  distribution  being  con- 
sidered) obtained  from  the  probability  of  damage  matrices,  this  section  of 
the  program  computes  the  probability  that  each  target  of  the  array  is 
damaged  for  each  volume  of  fire,  for  each  HOB.  The  probabilities  can  be 
averaged  over  a specified  HOB  distribution.  Using  the  average  probabilities 
(over  the  HOB  distribution)  that  each  target  sustain  damage,  the  probabil- 
ities that  various  numbers  of  targets  sustain  damage,  for  each  specified 
volley  size,  are  computed.  Also,  an  optimum  HUb  can  be  determined,  which 
is  then  treated  as  an  HOB  distribution  of  zero  width  to  obtain  the  quanti- 
ties discussed  above. 

Section  IV  identifies  those  fragments  of  the  fragment  distribution 
which  are  capable  of  perforating  a specified  plate  assuming  no  air  drag. 

For  each  such  fragment,  the  mass  and  velocity  of  the  most  massive  residual 
fragment  are  computed.  Fur  each  most  massive  residual  fragment,  the 
probability  of  incapacitating  a soldier,  given  that  he  is  hit,  is  computed 
using  the  methodology  of  reference  5.  The  average  probability  of  incapaci- 
tation over  the  distribution  of  perforating  fragments  is  computed. 

3.  MAJOR  ASSUMPTIONS 

a.  The  perforation  of  the  specified  armor  plates  by  the  specified 
fragments  is  described  by  the  methodology  of  Thor  reports  47  and  70  (ref- 
erences 1 and  2).  This  assumption  is  fairly  easy  to  change  should  a better 
methodology  become  available. 

b.  The  lightly  armored  structure  is  described  by  a set  of  flat, 
rectangular  armor  plates.  If  a plate  is  horizontal,  its  edges  must  be 
parallel  to  the  range  and  deflection  axes.  If  a plate  is  non-horizontal, 
it  can  deviate  from  a congruent,  but  horizontal,  plate  by  not  more  than  a 
rotation  downward  about  any  one  of  the  horizontal  plate's  edges. 

c.  The  probability  of  sustaining  at  least  N perforations  is  a 
good  measure  of  the  degree  to  which  indirect  artillery  fire  with  HE  pro- 
jectiles is  able  to  defeat  a lightly  armored  structure. 


4.  PROGRAM  SECTION  I - VULNERABILITY  OF  ARMOR  PLATES 
4. 1 The  Program. 

The  purpose  of  Section  I is  to  determine,  for  one  set  of  projectile 
fragmentation  data  and  one  set  of  projectile  coordinates,  the  expected 
number  of  perforations  sustained  by  one  armor  plate  of  given  material  and 
geometry  for  every  location  of  interest  of  that  plate  in  a horizontal  plane. 
This  information  is  placed  in  the  expected  perforations  matrix  file.  During 
the  computation,  various  outputs  are  generated,  which  are  of  interest  in  the 
themselves  and  which  may,  in  fact,  be  the  only  reason  for  processing  certain 
plates.  Conceptually,  the  methodology  of  this  section  is  straightforward. 

4.1.1  Horizontal  Plates.  Rather  than  compute  directly  the 
expected  number  of  perforations  sustained  by  the  target  plate  for  every 
location  of  interest,  we  consider  an  infinite  plate  of  the  specified  material 
and  thickness.  We  establish  a grid  of  square  cells  thereon,  and  fix  the 
grid's  location  with  respect  to  the  projectile  by  defining  its 
burst  point  (range,  deflection)  in  terms  of  grid  columns  and  rows.  The 
geometry  is  shown  in  Figure  4-1,  where: 

Z is  the  burst  point. 

OZ  is  the  vertical  separation  between  burst  point  and  plate. 

O'Z  is  the  HOB. 

ZQ  is  the  tangent  to  the  projectile  trajectory  at  the  burst  point. 

P is  the  point  currently  being  considered  for  fragmentation  damage, 
w is  the  angle  of  fal 1 . 

v is  the  supplement  of  the  polar  angle  to  the  point  P from  the  point  Z. 
n is  an  outward  drawn  unit  vector  normal  to  the  plate. 

C is  the  obliquity  angle  and  is  equal  to  'V . 

<p  is  the  azimuthal  angle  to  the  point  P. 

0 is  the  angle  between  the  projectile  axis  (assumed  to  coincide  with 
ZQ)  and  the  line  ZP. 


The  equation  for  the  angle  between  any  two  vectors  in  a spherical  polar 
coordinate  system  yields  equation  4-1  which,  in  a Cartesian  coordinate 
system,  reduces  to  equation  4-2,  which  is  used  to  determine  the  fragment 
spray  within  which  point  P lies. 


8 


cos  9 = cos  i|i  sin  <u  + sin  v cos  u cos  <}>  Equation  4-1 

cos  0 = (X^  + + Z^)  ^ (Z  sin  u + X cos  w)  Equation  4-2 

The  fragment  spray  within  which  the  point  P lies  is  determined,  and  the 
entire  cell  within  which  it  is  embedded  is  considered  to  be  covered  by  that 
spray.  Hence,  it  is  important  for  the  grid  cell  at  p)int  P to  subtend  a 
solid  angle  at  point  " which  is  small  with  respect  to  the  fragment  spray 
angle.  On  the  other  hand,  since  a computer  has  a limited  storage  capability, 
the  grid  must  be  limited  in  size.  Therefore,  the  grid  cell  should  be  large 
enough  to  produce  a quasi-infinite  plane  of  sufficient  size  to  encompass 
the  entire  area  over  which  significant  perforation  effects  may  occur. 

The  static  fragmentation  data  format  currently  used  breaks  the 
fragmentation  down  into  various  spray  angles  measured  from  the  projectile 
nose  and,  for  each  spray,  provides  a fragment  mass  distribution  and  the 
maximum  and  median  fragment  velocities  observed  within  that  spray. 

References  6 and  7 provide  examples  of  the  data  format  and  an  explanation 
of  the  data  collection  and  reduction  techniques.  The  fragment  sprays 
are  vectored  forward  to  account  for  the  projectile  velocity.  Because  the 
perforation  capability  of  a fragment  is  not  a linear  function  of  its  vel- 
ocity, and  also  because  consideration  of  each  individual  fragment  velocity 
(to  whatever  extent  such  data  are  available)  would  have  been  prohibitive  in 
terms  of  computer  time,  the  following  procedure  was  adopted  as  an  approxi- 
mation. The  number  of  perforations  sustained  by  the  (i,j)™  cell  is 
computed  assuming  that  all  fragments  have  the  appropriate  maximum  velocity 
and  then  recomputed  assuming  that  all  fragments  have  the  appropriate  median 
velocity.  The  expected  number  of  perforations  for  cell  (i,j)  is  approximated 
by  the  average  of  these  two  results.  Should  a comparison  of  these  approxi- 
mate values  with  experimental  data  indicate  significant  differences,  other 
combinations  of  the  two  computed  values  could  be  used  to  obtain  a better 
approximation. 

Once  the  fragmentation  data  have  been  vectored  and  the  spray(s) 
within  which  the  point  P lies  have  been  identified  (vectoring  the  sprays 
can  cause  them  to  overlap),  the  number  of  perforations  sustained  by  the 
cell  at  point  P is  computed  by  use  of  equation  4-3,  where 


M is  fragment  mass. 

V is  fragment  speed. 

Q is  either  mass  or  speed. 

r is  a subscript  identifying  residual  mass  or  residual  speed  of  the 
most  massive  fragment  particle  resulting  from  a perforation. 

M < 0 or  Vp  £ 0 means  no  perforation  occurs. 


t 


s is  a subscript  identifying  striking  mass  or  striking  speed  of  an 
impinging  fragment. 

-3/2 

K is  the  fragment  shape  factor  (reference  7);  K = m A where  A 
is  the  average  fragment  impact  area. 

t is  the  plate  thickness. 

C is  the  obliquity  angle  at  impact  of  the  fragment  with  the  plate. 

Cq,  cxq , Bq , Yg,  Xg  are  empirical  constants. 

Cn  an  -(2/3)aQ  (2/3)a«  + Br>  Yn  Xq 

Qr=Qs-10gt  gK  4 Ms  4 4 (sec  5)  4 V$  4 Equation  4-3 


This  equation,  from  Thor  Report  47  (reference  1).  yields  the  residual  mass 
and  speed  of  a fragment  after  perforation.  In  the  situation  at  hand,  that 
of  a horizontal  plate^cos  c is  just  cos  ¥ which  is  merely  the  ratio  of 
length  ZTF  to  length  ZP\ 

Thor  Report  70  (reference  2)  provides  equations  which  can  be 
used  to  extend  Equation  4-3  to  some  aluminum  alloys  not  considered  in 
Thor  Report  47  (reference  1).  These  equations,  after  some  algebraic 
manipulation  yield 


Qr(AL)  = i _ 


kQ(AL) 


kn(AL)  Q„(AL  2024) 


kg(AL’ 2024)  kg(AL  2024) 


Equation  4-4 


where  the  notation  is  the  same  as  in  equation  4-3  with  the  addition  of: 
kg  is  an  empirical  constant. 

AL  references  some  particular  aluminum  alloy. 


AL  2024  references  2024  aluminum. 

Reference  6 provides  equation  4-5  which  is  used  to  compute 
velocity  decay  due  to  air  drag: 

v$  = V0  exp  (-AA  • R • M$'(1/3))  Equation  4-5 


where 

Vg  is  the  velocity  with  which  the  fragment  strikes  the  plate 

* V is  the  velocity  with  which  the  fragment  is  projected  from  the 

fragment  source  point. 

4i 


n 


AA  is  the  air  drag  constant  provided  in  the  fragmentation  data. 

R is  the  distance  from  the  fragment  source  point  to  the  center  of 
the  plate  cell  being  struck. 

M is  the  mass  of  the  fragment  as  it  is  projected  from  the  fragment 
s source  point  and  as  it  strikes  the  plate. 

When  fragments  capable  of  perforating  the  plate  at  point  P (M  and 
Vr  both  greater  than  zero)  are  identified,  their  densities  per  unit  solid 
angle  are  summed  and  then  multiplied  by  the  solid  angle  subtended  by  the 
cell  at  point  P.  The  solid  angle,  n,  is  computed  using  equation  4-6, 
where  A is  the  surface  area  of  the  cell,  c is  the  obliquity  angle,  and  U 
must  be  smal 1 . 


P 

n = A • (X2  + Y2  + Z2)-1  • cos  ^ 


Equation  4-6 


The  preceding  is  done  for  both  the  maximum  and  median  velocities  of  the 
fragments  striking  point  P and  the  average  of  the  two  results  is  the 
expected  number  of  perforations  sustained  by  the  cell  at  point  P.  This 
procedure  is  performed  for  each  cell  in  the  grid. 

When  a direct  hit  or  very  nee*’  miss  occurs,  the  preceding 
computations  are  bypassed  and  a very  large  number  of  perforations  (cur- 
rently lCr)  is  inserted  into  the  matrix  elements  corresponding  to  the  quasi- 
infinite plate  cells  sustaining  a direct  hit  or  very  near  miss.  The  direct 
hit  problem,  then,  boils  down  to  identifying  the  cells  sustaining  a direct 
hit  or  very  near  miss.  This  is  done  in  the  following  manner. 

A sphere  with  radius  equal  to  the  blast  parameter  is  generated 
about  the  burst  point.  The  upper  hemisphere  is  replaced  by  a right 
circular  cylinder  whose  directrix  is  the  circle  formed  by  the  hemisphere 
edge.  The  intersection  of  the  resulting  solid  figure  with  the  quasi- 
infinite plate  is  a circle.  On  the  quasi-infinite  plate  is  constructed  the 
smallest  square  which  completely  encompasses  all  cells,  any  parts  of  which 
were  enclosed  by  the  circle.  The  cells  encompassed  by  this  square  are  the 
ones  considered  to  sustain  direct  hits  or  very  near  misses. 

The  expected  perforations  matrix  described  above  is  not  the 
desired  matrix;  it  merely  shows  the  expected  perforations  sustained  by  a 
horizontal,  square  target  of  unit  size  (one  cell)  for  every  position  of  the 
target.  To  obtain  the  expected  perforations  matrix  for  the  actual  target 
plate,  we  superimpose  on  the  actual  target  plate  a grid  having  the  same 
cell  size  and  cell  orientation  as  the  quasi-infinite  grid.  The  center  cell 
of  this  actual  target  plate  is  defined  by  equations  4-7: 

MIDX  = (NX©2)  + 1 

MIDY  = (NY  ©2)  + 1 


Equations  4-7 


where 

NX  is  the  number  of  cells  of  the  actual  target  plate  along  the  range 
axis. 

NY  is  the  number  of  cells  of  the  actual  target  plate  along  the 
deflection  axis- 

I 

©is  an  operation  defined  as  divide  and  then  truncate  the  fractional 
part. 


Defining  the  location  of  the  actual  target  plate  on  the  quasi-infinite 
plate  as  the  location  of  its  central  cell  on  the  quasi-infinite  plate, 
we  position  the  central  cell  of  the  actual  target  plate  successively  over 
each  cell  of  the  quasi-infinite  plate  and,  for  each  such  position,  we  sum 
the  expected  perforations  of  every  cell  covered  by  the  actual  target  plate. 
The  result  is  the  expected  perforations  matrix  for  the  actual  target  plate. 
This  matrix  is  now  multiplied  by  a scale  factor,  if  desired,  to  account  for 
any  slight  mismatch  in  surface  areas  that  may  exist  between  the  actual 
target  plate  and  the  plate  being  modeled. 

At  this  point,  the  expected  perforations  matrix  is  stored  in  the 
expected  perforations  matrix  file  and  various  types  of  information  may  be 
output,  as  will  be  discussed  in  paragraph  4.2. 

4.1.2  Non-Horizontal  Single  Plates.  The  plates  considered  here 
are  those  obtained  from  a congruent,  but  horizontal,  plate  by  a downward 
rotation  from  the  horizontal  about  only  one  of  the  horizontal  plate's  four 
edges.  Figure  4-2  shows  two  such  plates  rotated  about  the  deflection  axis. 
Figure  4-3  shows  two  such  plates  rotated  about  the  range  axis.  The 
following  sign  conventions  and  constraints  hold: 

(1)  0°  < fa | < 180° 

(2)  0°  <.  1 0 1 < 180° 

(3)  a is  positive  (negative)  when  a vector  drawn  outward  and 
normal  to  the  plate  has  a component  in  the  direction  of 
increasing  (decreasing)  range. 

(4)  B is  positive  (negative)  when  a vector  drawn  outward  and 
normal  to  the  plate  has  a conponent  in  the  direction  of 
increasing  (decreasing)  deflection. 

(5)  One  or  the  other  of  a and  B must  be  zero. 

In  thi'  program,  as  in  field  artillery  usage,  deflection  left  is  increasing 
and  deflection  right  is  decreasing.  The  treatment  of  such  plates  is 
similar  to  the  treatment  of  the  horizontal  plate,  but  there  is  an  increase 
in  complexity. 
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VIRTUAL  PLATES 


DEFLECTION 


RANGE 


Figure  4-2.  (U)  Single  Plates  With  a Rotation  from  the 
Horizontal  About  the  Deflection  Axis. 


The  target  plate  grid  no  longer  meshes  with  the  grid  of  the 
horizontal  plane;  therefore,  an  additional,  "virtual"  plate  is  appended 
to  the  top  edge  of  the  actual  target  plate.  Note  that  this  virtual  plate 
always  extends  in  the  direction  of  increasing  range  or  deflection  from  the 
top  edge  of  the  actual  plate.  The  central  cell  of  the  actual  target  plate 
is  defined  as  the  central  cell  of  the  virtual  plate,  and  it  is  determined 
by  applying  equations  4-7  to  the  virtual  plate. 

Whereas  in  paragraph  4.1.1,  we  considered  one  quasi -infinite 
plate  and  computed  the  expected  perforations  for  each  of  its  cells,  now  we 
consider  a series  of  quasi-infinite  virtual  plates  embedded  in  the  horizon- 
tal plane.  An  actual  target  strip  of  cells  (such  a strip  for  the  target  of 
Figure  4-2a  is  shown  in  Figure  4-4)  is  associated  with  each  of  the  virtual 
cells  for  a given  virtual  plate.  The  expected  number  of  perforations  for 
each  cell  of  this  strip  is  computed,  they  are  summed  to  obtain  the  total 
expected  perforations  for  the  strip,  and  this  quantity  is  stored  in  the 
matrix  as  the  expected  perforations  for  the  (i.j)th  virtual  cell.  It  must 
be  noted  that  computation  time  increases  rapidly  as  the  target  strip 
increases  in  length  for  any  given  cell  size.  The  equations  of  paragraph 
4.1.1  still  hold,  except  that  c is  no  longer  the  same  as  v.  Hence,  we  must 

us-  equation  4-8,  with  one  or  the  other  of  a and  6 always  zero. 

2 2 2 -("[/?)  Equation 

cos  S = (X  + Y + Z ) v '(Z  cos  a cos  3 + X sin  a + Y sin  3)  4-8 

Furthermore,  Z is  no  longer  constant  but  varies  as  we  sweep  down  the  strip 
of  cells. 

For  direct  hits  or  very  near  misses,  we  proceed  as  in  paragraph 
4.1.1  and  construct  a right  circular  cylinder  with  hemispherical  lower  end 
so  that  the  burst  point  is  at  the  center  of  the  parent  sphere.  The  virtual 

cells  whose  associated  actual  target  strips  are  enclosed  or  partially 

enclosed  by  the  cylinder-hemisphere  are  identified.  On  the  virtual  plate 
surface  the  smallest  rectangle  which  completely  encloses  all  such  virtual 
cells  is  generated.  The  virtual  cells  enclosed  by  this  rectangle  are 
identified,  and  their  associated  actual  target  strips  are  the  ones 
considered  to  sustain  direct  hits  or  very  near  misses. 

As  in  paragraph  4.1.1,  a sum  over  the  virtual  cells  is  performed 
to  obtain  the  expected  perforations  matrix  for  the  actual  target  plate. 

Of  course,  this  sum  is  now  1-dimensional  rather  than  2-dimensional. 

4.1.3  Non-Horizontal  Plate  Pairs.  The  concept  of  a non-horizon- 
tal  plate  pair  arises  from  an  examination  of  the  two  plates  shown  in  Figure 
4-2,  or  in  Figure  4-3,  and  from  the  realization  that  if  either  of  these 
plates  is  present  in  a lightly  armored  structure,  it  may  well  be  accompanied 
by  the  other  (e.g.,  the  two  sides  of  an  ARC).  For  this  case,  the  two  plates 
are  considered  as  one  plate  pair  and  the  symmetry  of  the  situation  is 
utilized  to  expedite  the  calculation.  Except  for  the  top  spacer  plate, 
nothing  new  is  introduced.  The  utility  of  the  concept  is  two-fold;  it  is 
faster  on  the  computer  than  doing  the  two  plates  separately,  and  it 
expedites  (for  the  machine)  and  facilitates  (for  the  user)  some  of  the 
subsequent  computations  of  program  section  II. 


Figure  4-5  shows  the  two  plates  of  Figure  4-2  combined  to  form  a 
plate  pair.  The  horizontal  plate  is  not  allowed  to  sustain  perforations 
and  serves  merely  as  a spacer.  The  central  cell  of  the  plate  pair  is 
defined  as  the  central  cell  of  the  top  spacer  and  is  computed  from  equations 
4-7  applied  to  the  top  spacer.  There  is  no  longer  a need  for  the  virtual 
plate  of  paragraph  4.1.2  unless  the  top  edges  of  the  two  non-horizontal 
plates  are  adjacent.  In  this  case,  the  geometry  can  be  visualized  by  letting 
the  virtual  plates  of  Figure  4-2  (or  Figure  4-3)  overlap  each  other  com- 
pletely. One  of  the  angles,  a and  8,  is  always  zero,  the  other  is  always 
positive  and  varies  between  0°  and  180°. 

4.1.4  Problems  Caused  by  Finite  Grid  Size.  In  order  to  obtain 
the  expected  perforations  matrix  for  an  actual  target  plate,  we  considered 
a quasi-infinite  plate  parallel  to  the  ground  plane  and  inscribed  with  a 
grid.  The  grid  is  currently  100  x 100  cells;  we  assumed  that  this  grid  was 
large  enough  so  that  the  resulting  expected  perforations  matrix  for  the 
actual  target  plate  would  cover  all  the  perforation  sustaining  locations 
of  the  actual  target  plate.  In  practice,  the  following  three  situations 
can  occur.  Although  the  discussion  addresses  only  horizontal  plates,  the 
generalization  to  other  plates  should  be  obvious. 

The  first  situation  is  that  tne  quasi-infinite  plate  is 
sufficiently  large  that  no  perforations  occur  beyond  its  limits.  Further- 
more, there  is  a border  (whose  required  width  will  become  apparent)  just 
within  each  edge  of  the  quasi -infinite  plate,  which  sustains  no  perforations. 
In  this  case,  there  is  no  problem.  The  summations  over  the  quasi-infinite 
plate  are  performed  without  errors  (henceforth  called  edge  effects)  due  to 
actual  target  cells  extending  beyond  the  edges  of  the  quasi-i nf i ni te  plate 
and  still  sustaining  perforations.  Furthermore,  by  considering  locations  of 
the  actual  target's  center  which  lie  beyond  the  limits  of  the  quasi -infi ni te 
plate,  we  see  that,  even  though  the  actual  target  extends  back  into  the 
quasi-infinite  plate,  it  extends  only  into  the  no  perforation  border.  This 
is  the  most  satisfactory  situation  because  the  final  result,  the  expected 
perforations  matrix,  contains  all  the  perforation  information  for  wery 
location  of  the  actual  target.  The  likelihood  of  occurrence  for  this 
situation  is  increased  as  the  cell  size  is  increased  and,  in  general, 
decreased  as  the  vertical  separation  between  the  burst  point  and  the  quasi- 
infinite plate  is  increased. 

The  second  situation  occurs  when  the  quasi-infinite  plate  is 
sufficiently  large  that  no  perforations  occur  beyond  its  limits,  but  the 
no  perforation  border  of  the  preceding  situation  either  is  absent  or  is  not 
wide  enough.  An  error  called  matrix  truncation  occurs  (there  are  locations 
not  represented  by  the  expected  perforations  matrix  for  which  the  actual 
target  will  sustain  perforations),  but  there  are  still  no  edge  effects. 

This  problem,  if  significant,  can  be  eliminated  by  solving  the  same  actual 
target  again  with  a different  projectile  location  with  respect  to  the  grid 
and  manually  piecing  together  the  several  expected  perforations  matrices, 
or  by  depending  upon  symmetry. 
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The  third  situation  occurs  when  the  quasi -infinite  plate  is  not 
sufficiently  large  to  preclude  perforations  beyond  its  boundaries.  Edge 
effects  occur  and,  of  course,  so  does  matrix  truncation. 

4.1.5  Symmetry.  The  program  assumes  that  the  fragmentation  data 
are  symmetric  across  the  range  axis  and  asymmetric  across  the  deflection 
axis.  Although  the  fragmentation  data  are  in  fact  symmetric  across  the  de- 
flection axis  for  certain  munitions  and  angles  of  fall,  this  symmetry  is  not 
utilized  by  the  program.  If  the  actual  target  plate  also  has  reflection 
symmetry  across  the  range  axis,  them  the  program  will  attempt  to  utilize  the 
reflection  symmetry  of  the  system  to  eliminate  any  edge  effects  which  may 
occur  across  the  low  deflection  edge  of  the  quasi-infinite  plate.  Thus,  in 
such  a case,  it  is  beneficial  to  locate  the  range  axis  (over  which  the  burst 
point  is  located)  at  the  low  deflection  edge  of  the  quasi-infinite  plate. 

The  following  example  illustrates  the  ideas  of  paragraphs  4.1.4  and  4.1.5. 

The  fragmentation  data  are  symmetric  across  the  range  axis  (as 
always)  and  across  the  deflection  axis  (angle  of  fall  = 90°).  The  actual 
target  plate  is  horizontal  and  so  enjoys  reflection  symmetry  across  both 
axes.  It  is  36  inches  long  (along  the  range  axis)  and  20  inches  wide  (along 
the  deflection  axis).  Considering  the  plate  thickness  and  material, 
vertical  separation  between  the  plate  and  the  burst  point,  and  the  fragmen- 
tation data,  we  expect  the  quasi-infinite  plate  to  sustain  perforations 
throughout  a circular  region  of  400  inch  radius  centered  under  the  burst 
point.  We  choose  a four-inch  cell  and  position  the  burst  point  at  Z as 
shown  in  Figure  4-6.  Because  of  symmetry,  there  are  no  edge  effects  across 
the  range  axis,  and  the  expected  perforations  matrix  can  be  reflected 
across  this  axis  to  eliminate  matrix  truncation.  Very  significant  edge 
effects  occur  in  columns  1 through  4;  however,  the  user  can  simply  delete 
these  columns  and  reflect  the  matrix  across  the  deflection  axis.  Matrix 
truncation,  but  not  edge  effects,  occur  across  the  high  deflection  edge  of 
the  matrix,  but  the  error  should  be  small  because  this  area  is  at  the  far 
edge  of  the  region  in  which  perforations  occur.  Roth  matrix  truncation  and 
edge  effects  occur  across  the  high  range  edge  of  the  matrix  but,  again,  the 
error  should  be  small. 

4.2  Output. 

a.  The  input  data  are  printed  out. 

b.  The  non-zero  portion  of  the  expected  perforations  matrix  is 
written  on  the  expected  perforations  file  with  a unique  record  number,  which 
is  the  ordinal  for  that  record.  These  records  are  sequentially  and  con- 
secutively numbered  for  all  cases  within  a batch  (paragraph  4.3).  Not  more 
than  2500  records  can  be  written  on  this  file.  If  the  entire  matrix  is  zero, 
then  a small  zero  matrix  is  written.  If  desired,  the  non-zero  portion  of 
the  matrix  can  be  printed.  If  the  entire  matrix  is  zero,  its  printing  is 
suppressed.  When  printing  occurs,  each  element  is  printed  in  an  16  format 
after  being  multiplied  by  a power  of  10  (input  data)  and  truncated. 
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c.  The  non-zero  portion  of  the  probability  of  damage  matrix  for 
each  criterion  is  printed  out  if  desired.  If  the  expected  perforations 
matrix  is  zero,  none  of  its  probability  of  damage  matrices  are  computed. 

If  the  probability  of  damage  matrix  associated  with  the  expected  perfor- 
ations matrix  for  some  damage  criterion  is  zero,  its  printing  is  suppressed 
and  computation  of  any  succeeding  probability  of  damage  matrices  which  may 
have  been  requested  for  that  expected  perforations  matrix  is  suppressed. 

Each  element,  P.(i,j),  of  the  probability  of  damage  matrix  is  computed  from 
the  corresponding  element  of  the  expected  perforations  matrix  using  the 
specified  damage  criterion  and  assuming  a Poisson  distribution. 

d.  Lethal  area  (reference  3),  average  probability  of  damage  for 
all  locations  of  the  actual  target  plate  or  plate  pair  for  which  the 
probability  of  damage  is  non-zero,  and  the  ground  surface  area  over  which 
there  is  a non-zero  probability  of  damage  (henceforth  called  "damage  area"), 
are  printed  for  each  probability  of  damage  matrix  computed.  These  values 

do  not  include  multiplication  by  two  to  account  for  reflections  of  a half- 
infinite ground  plane  across  the  range  axis. 

e.  A contour  plot  of  each  probability  of  damage  matrix  computed 
is  printed  out  if  desired.  The  contour  plot  is  suppressed  if  the  expected 
perforations  matrix  or  the  probability  of  damage  matrix  is  null.  An  example 
is  at  Figure  4-7. 

f.  Whenever  a contour  plot  is  generated,  the  program  also 
generates  and  stores  a damage  function  on  the  damage  function  file  with  a 
unique  record  number,  which  is  the  ordinal  for  that  record.  The  records 
are  sequentially  and  consecutively  numbered  for  all  functions  generated 
within  a batch  (paragraph  4.3).  Not  more  than  5000  such  functions  can  be 
stored  on  this  file.  The  damage  function  generated  from  the  probability  of 
damage  matrix  of  Figure  4-7  is  at  Figure  4-8.  A region  is  that  portion  of 
the  matrix  for  which  the  P<j(i,j)  are  less  than  or  equal  to  one  contour  value 
but  greater  than  the  next  lower  contour  value.  A region  need  not  be  contin- 
uous. Rectangularized  regions  are  those  rectangular  regions  defined  by  the 
set  of  smallest  rectangles  completely  enclosing  the  set  of  actual  regions 

on  a one-for-one  basis.  The  statement  that  these  rectangularized  regions 
are  unnormalized  means  that  lethal  area  and  damage  area  are  not  conserved 
in  their  generation;  hence,  before  they  can  be  used,  these  damage  functions 
must  be  normalized.  The  various  regions  in  one  damage  function  are  written 
on  the  file  from  outermost  to  innermost.  The  maximum  number  of  such  regions 
is  11,  corresponding  to  the  regions  within  contour  lines  for  P . = 1 . , .9, 

.8,  . . . , .1 , and  .0005.  0 

4.3  Input  Data  Deck 

This  deck  inputs  a "batch",  or  "logical  grouping",  of  cases. 

All  cases  within  a batch  have  the  same  cell  size,  ground  plane  grid  size, 
projectile  location  with  respect  to  the  grid,  and  fragmentation  data.  For 
the  purposes  of  this  paragraph,  a "target"  is  one  plate  or  one  plate  pair. 

A "case",  then,  is  one  target  at  one  plate  thickness  and  one  projectile 
HOB.  Each  target  is  solved  for  every  specified  thickness  at  every 
specified  HOB. 
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Unless  otherwise  stated,  each  data  field  is  10  columns  wide  and 
each  numeric  datum  is  in  a real  format  with  an  explicit  decimal  point. 

Data  fields  not  specifically  addressed  are  ignored  and  may  be  used  for  any 
desired  purpose,  such  as  user  comments  within  the  data  deck. 

4.3.1  The  Data  Deck.  The  data  deck  for  program  section  I is  at 
Figure  4-9. 

4.3.2  Definition  of  Each  Input  Variable 


Card  1: 

Columns  1-10  - the  words  "BASIC  DATA". 

FILES  - the  unique  whole  number  identifying  the  first  perforation  matrix  to 
be  generated  in  this  batch  of  cases.  The  perforation  matrix  file 
cannot  contain  more  than  2500  records  (matrices).  FILES  can  be  any 
whole  number  from  1 thru  2500.  FILES  is  the  first  value  of  a con- 
secutive series  of  numeric  indices  for  the  matrices  to  be  generated 
and  filed  in  the  current  batch  of  cases.  FILES  itself,  however, 
need  not  be  consecutive  or  even  sequential  with  respect  to  the 
matrices  stored  in  the  perforation  file  by  all  preceeding  batches. 
If  FILES  is  not  greater  than  the  index  of  the  highest  indexed  matrix 
placed  in  the  file  by  all  preceeding  batches,  some  or  all  of  che 
earlier  matrices  may  be  overlaid  by  matrices  of  the  current  batch. 
There  is  a default  value  for  FILES.  If  FILES  is  blank  or  0.  its 
value  is  computed  to  be  one  greater  than  the  index  of  the  highest 
indexed  matrix  already  in  the  file.  If  the  file  is  empty,  the 
defaul t fcr  FILES  is  1 . 

CELL1  - the  edge  length,  in  inches,  of  the  square  cell  which  will  be  used 
to  grid  both  the  target  and  the  ground  plane. 

NXX,NYY  - the  numbers  of  cells  in  the  ground  plane  grid  (and  hence  in  the 

expected  perforations  matrix)  along  the  range  and  deflection  axes, 
respecti vely ; 1.  £ NXX  _<  100.;  1.  <_  NYY  £ 100.;  in  practice,  these 
variables  would  be  less  than  100.  only  to  save  machine  time. 

IC0L,IR0W  - related  to  the  range  ana  deflection  coordinates  of  the  burst 
point  in  terms  of  ground  plane  grid  columns  and  rows,  respect- 
ively; grid  columns  increase  with  range  and  grid  rows  increase 
with  deflection  left;  (IC0L,IR0W)  = (50., 0.)  would  yield  burst 
point  ground  plane  grid  coordinates  = (50. 5, .5),  corresponding 
to  a location  between  the  50th  and  51st  grid  columns  and  between 
the  zeroeth  (non-existent)  and  first  grid  rows;  in  practice,  the 
usefulness  of  these  variables  is  to  shift  the  ground  plane  grid 
location  with  respect  to  the  (fixed)  burst  point;  a burst  point 
outside  the  ground  plane  grid  is  allowed,  with  the  exception 
stated  in  paragraph  5.3. 
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Card  2: 


NHOB  - the  number  of  HOB's  to  be  considered  in  this  batch  of  cases; 

1.  < NHOB  < 25. 

HHHOB(I)-  the  Ith  HOB  to  be  considered  in  this  batch  of  cases;  in  meters. 

As  many  as  three  additional  cards  may  occur  here  to  input  up  to  25  values 
of  HHHOB;  each  additional  card  starts  with  a value  of  HHHOB  in  the  first 
data  field  (columns  1-10). 

Card  3: 

NTHCK  - the  number  of  armor  plate  thicknesses  to  be  considered  in  this 
batch  of  cases;  each  thickness  is  considered  for  each  HOB; 

1.  < NTHCK  < 25. 

i.  L 

TTHCK(I)  - the  Itn  thickness  to  be  considered  in  this  batch  of  cases;  in 
inches . 

As  many  as  three  additional  cards  may  occur  here  to  input  up  to  25  values 
of  TTHCK;  each  additional  card  starts  with  a value  of  TTHCK  in  the  first 
data  field  (columns  1-10). 


Card  4: 


NTRGT  - the  number  of  targets  to  be  considered  in  this  batch  of  cases; 

1.  £ NTRGT  _<  25.;  each  target  is  considered  for  each  thickness, 
thus,  the  number  of  cases  in  this  batch  is  NTRGT*NTHCK*NHOB;  an 
overall  constraint  is  l.<  NTRGT*NTHCK*NHOB  <2500. 

Card  5: 

There  are  NTRGT  cards  5;  each  provides  information  for  one  target  in  this 
batch  of  targets;  all  targets  will  be  solved  for  the  first  thickness  at  the 
first  HOB,  then  for  the  second  thickness  at  the  first  HOB,  and  so  on  until 
all  thicknesses  have  been  exhausted,  then  for  the  first  thickness  at  the 
second  HOB,  and  so  on  until  all  HOB's  have  been  exhausted. 

NXTGT(I)  - the  number  of  cel  1 s.^along  the  range  axis  in  the  I ^ plate  or  in 
the  spacer  of  the  plate  pair. 

NYTGT(I)  - the  number  of  cells  along  the  deflection  axis  in  the  I plate  or 
in  the  spacer  of  the  I plate  pair. 

• NZTGT(I)  - the  number  of  cells  along  the  slope  in  the  Ith  plate  or  plate 

pair;  zero  or  blank  for  horizontal  plates. 

. NSIDE(I)  - zero  or  blank  for  horizontal  plates;  1.  for  single  non-horizontal 

plates;  2.  for  plate  pairs. 

* 


ALPHA(I) 

BETA(I) 

TT  ( I ) 
SC(I) 

SH(I) 

Card  6: 
NKRIT 

KRIT(I) 

Card  7: 
Columns  1 

Card  8: 


- 0.  £ | ALPHA( I ) | < 180.;  the  angle  in  degrees  from  the  horizontal 
down  to  the  outer  surface  of  the  plate  or  plate  pair; 
positive  for  plate  pairs  and  plates  whose  outward  drawn  normal 
vector  has  a component  in  the  direction  of  increasing  range; 
negative  otherwise;  zero  for  plates  rotated  about  the  range  axis 
and  for  horizontal  plates. 

- 0.  j<  | BETA ( I ) | < 180.;  the  angle  in  degrees  from  the  horizontal 
down  to  the  outer  surface  of  the  Ith  plate  or  plate  pair; 
positive  for  plate  pairs  and  plates  whose  outward  drawn  normal 
vector  has  a component  in  the  direction  of  increasing  deflection 
(deflection  left);  negative  otherwise;  zero  for  plates  rotated 
about  the  deflection  axis  and  for  horizontal  plates. 

- the  height  of  the  top  surface  or  top  edge  of  the  I**1  plate  or 
plate  pair  above  the  ground  plane;  in  inches. 

- 5 column  field;  actually  SCALE ( I ) ; a scale  factor  used  to  adjust 
for  minor  mismatch  between  actual  and  model  plate  surface  areas; 
has  no  effect  on  direct  hit  contributions;  equal  to  the  surface 
area  of  the  plate  being  input  to  the  program  divided  by  the  sur- 
face area  of  the  plate  being  modeled;  default  value  is  1. 

- 5 column  field;  actually  SHOW(I);  non-negative  whole  number;  the 
power  to  which  10.  will  be  raised  and  then  used  to  multiply  each 
element  of  the  expected  perforations  matrix  prior  to  output  of 
that  matrix;  if  SH(I)  = 0.  or  blank,  the  Ith  matrix  is  not  output; 
the  output  field  for  each  element  of  the  matrix  is  16,  including 
spaces  between  elements. 


- the  number  of  damage  criteria  which  will  be  considered  for  each 
case  in  this  batch;  1.  £ NKRIT  £ 7. 

- the  Ith  damage  criterion;  KRIT(I)  > 0;  the  input  vector  must  be 
monotone  increasing. 


-10  - the  words  THOR  DATA;  identifies  the  following  two  cards  as 
providing  Thor  constants;  from  and  in  the  same  units  as  Thor 
reports  47,  70  (references  1,2);  if  not  present,  cards  8 and  9 
also  must  be  not  present  and  Thor  constants  from  the  immediately 
preceding  batch  of  cases  will  be  used. 


- 10  alphanumeric  characters ; the  user's  own  name  for  the  material 
to  which  the  Thor  constants  apply. 
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CV,  ALFAV,  BETAV,  GAMAV,  LAMDAV,  KFCTRlf  - Thor  constants  fof  €he  residual 
velocity  equations  with  no  particular  fragment  shape  assumed; 
KFCTRV  is  blank  or  0.  if  the  armor  is  not  an  aluminum  for  which 
K- factors  are  available,  as  in  Thor  report  70  (reference  2). 


Card  9: 

/ 

CM,  ALFAM,  BETAM,  GAMAM,  LAMDAM,  KFCTRM  - Thor  constants  for  the  residual 

mass  equations  with  no  particular  fragment  shape  assumed;  KFCTRM 
is  blank  or  0 . if  the  armor  is  not  an  aluminum  for  which  K- 
factcrs  are  available,  as  in  Thor  report  70  (reference  2). 

Card  10: 

W - the  projectile  angle  of  fall  in  degrees. 

VTERM  - The  projectile  terminal  velocity  in  feet  per  second. 


| 


i 
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BLAST  - radius  in  inches  of  the  blast  effects  sphere,  within  which  the 
probability  of  damage  is  1. 

Card  11: 

Columns  1-10  - the  word,  FRAGMENTAT ; identifies  subsequent  cards  as  providing 
the  fragmentation  distribution;  if  not  present,  cards  12  through 
17  also  must  not  be  present  and  the  fragmentation  data  from  the 
immediately  preceding  batch  of  cases,  or  from  the  immediately 
preceding  execution  of  program  Section  IV,  will  be  used. 

SHAPE  - the  fragment  shape  factor;  in  grains/cubic  inch. 

AA  - a constant  in  the  velocity  decay  equation  (see  equation  4-5). 
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Card  12: 

NSPRA  - the  number  of  fragment  sprays  in  the  fragmentation  data; 

3.  < NSPRA  < 46. 

NSP(I)  - the  number  of  mass  values  in  the  Ith  spray;  this  vector  may  be 
continued  on  as  many  as  5 additional  cards,  starting  in  column  1 
of  each,  for  NSPRA  elements  in  NSP;  1.  <_  NSP(I)  <_  40.,  all  I. 

Card  13: 

There  may  be  as  many  as  46  sets  of  as  many  as  11  cards  each  to  input  as  many 
as  40  Q(J,K),  M(J,K)  pairs  for  each  of  the  NSPRA  fragment  sprays;  K runs 
from  1 to  NSPRA;  J runs  from  1 to  NSP(K). 

K - the  spray  number. 


B 
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Q(J,K) 

M(J,K) 

Card  14: 

K 

SPRBS(K) 

VSTMN(K) 

VSTMX(K) 

NUM(K)  - 
SGST(K) 
Card  15: 

Col umns 

NSTAPE 


the  fraction  of  all  fragments  in  the  spray  residing  in  the 
jth 

mass  interval. 

the  mass .assigned  to  fragments  residing  in  the  mass  interval 
of  the  Ktn  fragment  spray;  in  grains;  for  any  given  K,  the  M(J,K) 
must  be  monotone  increasing. 


There  may  be  as  many  as  46  cards  #14  to  define  NSPRA  fragment 
sprays. 

- the  spray  number. 

, SPRBS(K+1J  - the  leading  and  trailing  edges  of  the  Kth  fragment 

spray;  in  degrees;  the  leading  (trailing)  edge  of  the 
first  (last)  spray  must  be  0.°(180.°);  the  sprays  must 
be  continuous  from  0°  through  180°. 

1.  L 

- the  median  fragment  velocity  at  the  midpoint  of  the  Ktn  fragment 
spray;  in  feet  per  second;  the  midpoints  of  the  first  and  last 
sprays  are  defined  as  0.°  and  180.°,  respectively. 

- the  maximum  fragment  velocity  at  the  midpoint  of  the  Kth  fragment 
spray;  in  feet  per  second;  the  midpoints  of  the  first  and  last 
sprays  are  defined  as  0.°  and  180.°,  respectively. 

iL 

the  number  of  fragments  in  the  Ku  fragment  spray;  if  0.  or  blank, 

will  be  computed  from  SGST(K). 

- the  density  of  fragments  (number  per  steradian)  in  the  Kth  fragment 
spray;  if  0.  or  blank,  will  be  computed  from  NUM(K). 


An  optional  card  which  may  use  either  of  two  formats. 

The  first  format  is: 

-10  - the  word  CONTOUR;  causes  a contour  map  of  the  probability  of 
of  damage  matrix  to  be  output  for  each  case  and  damage  criterion 
of  the  current  batch;  also  causes  generation  of  a damage  function, 
storage  thereof  on  the  damage  function  file,  and  printed  output 
thereof. 

the  unique  whole  number  identifying  the  first  damage  function  to 
be  generated  in  this  batch  of  cases.  The  damage  function  file 
cannot  contain  more  than  5000  records  (functions).  NSTAPE  can  be 
any  whole  number  from  1.  thru  5000.  NSTAPE  is  the  first  value  of 
a consecutive  series  of  numeric  indices  for  the  functions  to  be 
generated  and  filed  in  the  current  batch  of  cases.  NSTAPE  itself, 


however,  need  not  be  consecutive  or  even  sequential  with  respect  to 
the  functions  stored  in  the  damage  function  file  by  all  preceeding 
batches.  If  NSTAPE  is  not  greater  than  the  index  of  the  highest 
indexed  function  placed  in  the  file  by  all  preceeding  batches, 
some  or  all  of  the  earlier  functions  may  be  overlaid  by  functions 
of  the  current  batch.  There  is  a default  value  for  NSTAPE.  If 
NSTAPE  is  blank  or  0.,  its  value  is  computed  to  be  one  greater 
than  the  index  of  the  highest  indexed  function  already  in  the  file. 
If  the  file  is  empty,  the  default  for  NSTAPE  is  1. 


Columns  21-30  - the  words  NO  PROBABI;  causes  suppression  of  the  printing  of 
all  the  probability  of  damage  matrices  for  this  batch. 

The  alternate  format  is: 


Columns  1-10  - the  words  NO  PROBABI,  as  above. 


Columns  31-40  - the  word  CONTOUR,  as  above. 
Columns  41-50  - the  variable  NSTAPE,  as  above. 


31 
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5.  PROGRAM  SECTION  II  - VULNERABILITY  OF  ARMORED  STRUCTURE 


5.1  The  Program. 

The  function  of  program  Section  II  is  to  assemble  the  expected 
perforations  matrices  (generated  by  Section  I and  retrievable  from  the 
perforation  matrix  file)  for  as  many  as  seven  plates  or  plate  pairs  to  obtain 
the  expected  perforations  matrix,  probability  of  damage  matrices,  and  damage 
functions  for  one  complete  target.  All  the  matrices  which  are  combined  for 
any  one  case  must  be  for  the  same  cell  size  and  projectile  HOB,  angle  of 
fall,  and  terminal  velocity.  There  can  be  up  to  25  cases  in  one  batch. 

The  manner  in  which  the  expected  perforations  matrices  associated  with 
several  plates  are  combined  to  form  the  expected  perforations  matrix  for  a 
complete  target  is  best  described  through  the  use  of  an  example. 

5.1.1  Assembling  a Complete  Target.  Consider  the  target  shown  in 
Figure  5-1,  which  shows  a side  view  with  elevation,  a side  view,  and  a front 
view.  Plates  seen  on  edge  have  tick  marks  to  indicate  cell  boundaries. 

The  extent  of  plates  seen  on  edge  is  indicated.  The  scale  and  coordinate 
system  are  common  to  all  three  views.  Data  for  each  plate  or  plate  pair  are 
shown  in  Table  5-1. 

Plate  A is  embedded  in  the  spacer  BB  which  can  sustain  no  perfor- 
ations. Between  the  edges  of  plate  A and  the  top  edges  of  sides  B is  a 
region  which  could  correspond  to  heavily  armored  parts  of  an  actual 
structure.  The  shaded  region  is  the  virtual  plate  associated  with  plate  C. 
Reflection  symmetry  exists  across  the  range  axis. 

Figures  5-2,  5-3,  and  5-4  provide  examples  of  additional  plates 
which  could  be  part  of  the  complete  structure,  but  which  are  not  considered 
in  the  current  example.  Partial  data  are  given  in  Table  5-1.  These  plates 
are  shown  to  indicate  how  reasonably  complex  structures  can  be  built  up  and 
to  show  various  possible  values  of  ALPHA,  BETA,  AND  TT.  The  points  Q,  R, 
and  S of  Figures  5-2,  5-3,  and  5-4  correspond  to  those  same  points  in 
Figure  5-1. 

Consider  plate  A with  central  cell  D and  plate  pair  B with  central 
cell  E.  We  define  plate  A as  the  reference  plate  for  the  structure;  its 
central  cell  will  become  the  central  cell  for  the  entire  structure. 

Using  the  ground  plane  grid  of  Figure  4-6,  element  (I,J)  of  the 
appropriate  expected  perforations  matrix  is  the  expected  number  of  per- 
forations sustained  by  plate  A when  its  central  cell  (D)  is  positioned  over 
cell  (I,J)  of  the  ground  plane.  Similarly,  element  (I+1,J)  of  another 
expected  perforations  matrix  is  the  expected  number  of  perforations  sustained 
by  plate  pair  B when  its  central  cell  (E)  is  positioned  over  cell  (I+1,J)  of 
the  ground  plane.  Thus,  for  all  rows  J and  all  columns  I: 


Figure  5.1  (U)  Model  of  a Complete  Target. 
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Figure  5.4.  (U)  Rear  Plate  Which  Could  be  Added 
to  the  Model  of  figure  5.1. 


EPMAjB(I,J)  = EPMa(I,J)  + EPMb(I+1,J) 


Equation  5-1 


where  EPM^  B ( I , J ) means  element  (I,J)  of  the  expected  perforations  matrix 

for  the  structure  made  up  of  plate  A and  plate  pair  B. 

Similarly, 


EPMA,B,C(I,J)  = EPMA,B(I,J)  + EPMC(I+8’J) 


Equation  5-2 


EPMa  b c(I,J)  = EPMa(I,J)  + EPMB(I+1,J)  + EPMC(I+8,J)  Equation  5-3 


Thus,  the  expected  perforations  matrix  for  a complete  structure  is  obtained 
by  adding,  with  appropriate  offsets,  the  expected  perforations  matrices  for 
individual  components. 

It  would  be  most  efficient  to  first  compute  the  offsets  in  range 
and  deflection  and  then  add  all  the  matrices  together,  element  by  element, 
at  once,  as  indicated  by  Equation  5-3,  This,  however,  was  beyond  the  stor- 
age capability  of  the  computer  for  which  this  program  was  originally 
developed  (BRLESC  I / 1 1 of  the  USABRL);  therefore,  Equations  5-1  and  5-2  are 
used.  The  first  two  components  are  combined  to  obtain  an  intermediate 
structure  which  is  then  treated  as  a component  with  which  the  next  component 
is  combined,  and  so  on.  Therefore,  the  first  component  (the  one  identified 
by  logical  record  number  FILE(l) ) is  always  used  as  the  reference  component 
and  its  central  cell  always  becomes  the  central  cell  of  the  complete  struc- 
ture. For  reasons  which  will  become  apparent,  the  choice  of  the  order  in 
which  components  will  be  combined  is  a matter  of  some  importance. 

Only  those  matrix  addition  offsets  due  to  the  offsets  of  central 
cells  from  the  reference  central  cell  have  been  discussed.  There  are  two 
other  offsets  which  are  of  interest,  and  over  which  the  user  has  only 
partial  control. 

Consider  the  expected  perforations  matrices  associated  with 
plates  A and  C,  and  assume  (for  simplicity)  that  the  dynamic  fragmentation 
projected  from  the  burst  point  is  spherically  symmetric.  Since  maximum 
damage  occurs  when  the  struck  plate  is  perpendicular  to  the  trajectories 
of  the  striking  fragments  (obliquity  angle  = 0°),  it  is  clear  that  the 
region  of  maximum  damage  for  plate  C will  be  in  lower  numbered  columns  than 
it  will  be  for  plate  A,  for  a wide  range  of  HOB's.  Hence,  the  user  well 
may  have  positioned  the  projectile  over  different  grid  points  when  com- 
puting EPM.  and  EPMf  in  order  to  minimize  matrix  truncation  (paragraph 
4.1.1)  for  each  of  Lthese  plates.  The  program  automatically  compensates 
for  any  such  offsets  in  projectile  locations. 

The  following  describes  the  second  type  of  offset.  Assume 
NXX  = NYY  = 100  for  each  expected  perforations  matrix.  The  various  offsets 
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applied  when  combining  these  component  matrices  may  result  in  an  expected 
perforations  matrix  for  the  complete  structure  which  extends  beyond  one  or 
more  sides  of  the  NXX  x NYY  cell  grid.  Each  time  a component  matrix  is 
added  into  the  result  of  previous  additions  (or  into  the  reference  matrix, 
if  th^  addition  is  the  first  one  to  be  performed),  the  program  will  attempt 
to  minimize  matrix  truncation  by  rigidly  translating  the  entire  plates/ 
projectile  system  within  the  grid.  However,  the  program  always  gives 
preference  to  the  prior  matrices  over  the  later  matrices;  that  is,  regard- 
less of  the  magnitudes  of  the  matrix  elements  being  lost,  and  regardless  of 
the  number  of  such  elements,  the  program  will  never  shift  any  part  of  a 
prior  matrix  or  matrix  sum  out  of  the  grid  in  order  to  accommodate  the 
matrix  currently  being  added  in.  Therefore,  consideration  must  be  given 
to  the  order  in  which  the  component  matrices  are  presented  to  the  program 
(i.e.:  the  order  of  the  NFILE  elements  of  the  set  { FI LE ( I ) } ; see  paragraph 
5.4)  and  the  dominant  matrices  should  be  presented  first. 

5.1.2  Symmetry  and  an  Example.  As  in  the  case  for  Section  I of 
the  program,  Section  II  can  make  use  of  symmetry  across  the  range  axis,  and 
the  user  can  utilize  symmetry  across  the  deflection  axis.  Consider  the 
structure  of  Figure  5-1  as  a target  oriented  in  the  direction  of  increasing 
range.  Clearly,  reflection  symmetry  exists  across  the  range  axis  for  each 
target  component  and  also  for  the  complete  target.  Therefore,  the  burst 
point  for  each  component  should  be  located  between  ground  plane  grid' rows  0 
(non-existent)  and  1 so  that  each  component's  expected  perforations  matrix 
will  contain  maximum  information.  As  the  target  is  assembled,  there  will 
be  no  need  to  shift  the  burst  point  from  this  zero  deflection  location 
(shifts  in  range  may  occur)  and  the  final  matrix  also  will  contain  maximum 
information.  A similar  discussion  applies  to  the  same  vehicle  oriented  in 
the  direction  of  decreasing  range. 

Rotate  the  structure  of  Figure  5-1  clockwise  by  90°  and  consider 
it  as  a target  oriented  across  the  range  axis  in  the  direction  of  decreasing 
deflection  (Figure  5-5).  (Note  that  the  virtual  plate  belonging  to  plate  C 
has  changed  its  orientation  with  respect  to  plate  C,  and  now  lies  on  top  of 
one  end  of  plate  A.)  Since  the  fragmentation  is  symmetric  across  the  range 
axis,  it  is  clear  that  the  solution  over  the  entire  ground  plane  differs 
from  that  for  the  same  target  oriented  in  the  direction  of  increasing 
deflection  only  by  a reflection  across  the  range  axis.  Thus,  we  need  only 
one  of  the  two  orientations  for  the  target  broadside  to  the  range  axis.  We 
cannot,  however,  proceed  (as  in  the  preceding  paragraph)  to  solverover  the 
half-infinite  ground  plane,  and  depend  on  reflection  symmetry  of  the  target 
itself  to  obtain  the  complete  solution,  because  target  reflection  symmetry 
does  not  exist.  The  problem  can  be  overcome  in  either  of  two  ways.  The 
first  is  to  position  the  burst  point  somewhere  over  the  center  of  the  ground 
plane  grid.  This  is  completely  satisfactory  so  long  as  the  user  can  stand 
the  (extreme)  matrix  truncation  which  may  occur.  The  other  method  will  be 
discussed  in  paragraph  5.3.  At  this  point,  we  present  an  example  to  illus- 
trate the  solution  over  the  half-infinite  plane  of  the  target  facing  in  the 
direction  of  decreasing  deflection. 


We  have  stored  in  a file  the  expected  perforations  matrices  for 
each  target  component  (plate  A,  plate  pair  B,  and  plate  C of  Figure  5-5). 
When  these  matrices  were  computed,  we  chose  the  projectile  burst  locations 
with  some  care;  for  components  A and  B,  it  was  between  grid  rows  0 (non- 
existent) and  1,  and  for  component  C,  it  was  between  grid  rows  9 and  10. 
Plate  A is  considered  the  reference  plate  and  its  central  cell  is  the 
central  cell  for  the  entire  structure.  As  stated  previously,  the  perfor- 
ation matrices  for  the  several  components  are  additive  so  long  as  the  proper 
offsets  are  maintained.  To  determine  which  element  of  a component's 
perforation  matrix  should  be  added  to  an  element  of  the  reference  plate's 
or  reference  structure's  perforation  matrix,  we  use  Equation  5-4  and  a 
similar  equation  for  matrix  columns. 

(component  matrix  row  to  be  added  in) 


= (structure  matrix  row  under  consideration) 

+ (row  offset  of  component  central  cell  with  respect  to  structure 
central  cell) 


+ (row  offset  of  burst  point  in  component  matrix  with  respect  to 


burst  point  in  structure  matrix) 


Equation  5-4 


For  plate  pair  B,  the  last  two  terms  of  Equation  5-4  are  -1  and  0,  respect- 
ively. Furthermore,  there  are  no  column  offsets.  Therefore, 


EPMA  b(I,J)  = EPMa(I,J)  + EPMb(I,J-1)  Equation  5-5 

We  note  that  when  cell  D is  in  grid  row  1,  there  is  no  information  available 
for  plate  pair  B and  hence  we  would  have  an  edge  effect  error,  except  for 
the  fact  that  plate  pair  B enjoys  reflection  symmetry  across  the  range  axis. 


A general  statement  is  that,  if  matrix  B has  central  cell  row 
offset,  K,  and  burst  point  row  offset,  L,  then 

EPMa  b(I,J)  = EPMa(I,J)  + EPMb(I,J+K+L)  Equation  5-6 

Furthermore,  if  component  B enjoys  reflection  symmetry  across  the  range  axis 
and  if  K < 0 and  K + L < 0,  then,  when  J £ |K  + L| , Equation  5-6  is  replaced 
by  Equation  5-7. 


EPMa  b(I,J)  = EPMa(I,J)  + EPMb(I,1-(J+K)+L)  Equation  5-7 

Returning  to  the  specific  case  at  hand,  we  have  Equation  5-5  which  holds 
except  when  cell  D is  in  grid  row  1.  When  cell  D is  in  grid  row  1,  we  use 
Equation  5-8,  obtained  from  Equation  5-7. 

EPMa  b(I,J)  = EPMa(I,J)  + EPMg(1 ,2-J ) Equation  5-8 
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Thus  far,  we  have  obtained  the  expected  perforations  matrix  for  the 
structure  of  Figure  5-5,  not  including  plate  C.  The  structure's  central 
cell  is  cell  D.  The  range  axis  is  still  located  between  grid  rows  0 (non- 
existent) and  1.  The  structure  does  not  enjoy  reflection  symmetry  across  the 
range  axis  (because  plate  pair  B sticks  out  in  front  of  plate  A,  but  not  in 
back  of  it),  but  this  is  of  no  importance  because  the  structure  is  not  added 
into  anything;  rather,  components  are  added  into  the  structure.  To  add 
plate  C into  the  structure,  we  use  equation  5-9,  obtained  from  equation  5-6 
with  K = -9  and  L = 9 


Vb,C(I'J>  ■ epma,b(i'j>  + EPMc(I-J> 


Equation  5-9 


It  should  now  be  clear  that  the  entire  purpose  for  locating  the 
range  axis  for  plate  C between  grid  rows  9 and  10  was  to  avoid  edge  effects 
due  to  plate  C when  the  central  cell  of  the  structure  is  near  the  burst 
point.  Equation  5-7  could  not  have  been  used  because  plate  C does  not 
enjoy  reflection  symmetry  across  the  range  axis. 

5.2  Output. 

Program  status  at  this  point  is  very  similar  to  its  status  after 
computations  for  a target  component  have  been  completed  by  Section  I.  The 
output  is  the  same  as  in  paragraph  4.2  with  the  following  exceptions: 

a.  The  expected  perforations  matrix  for  the  complete  structure 
is  not  stored  on  the  perforation  matrix  file. 

b.  There  is  no  suppression  of  probability  of  damage  matrices  or 
of  the  contour  plots  and  damage  functions  due  to  null  matrices. 

5 . 3 An  Additional  Capability  for  Targets  Displaying  Asymmetry  Across 

the  Range  Axis. 

Consider  the  target  of  Figure  5-5  and  allow  it  to  assume  both 
orientations  perpendicular  to  the  range  axis  in  both  half-infinite  parts  of 
the  ground  plane.  Figure  5-6  illustrates  the  four  possibilities.  For  each 
case,  the  target  is  symmetric  across  the  deflection  axis  and  asymmetric 
across  the  range  axis.  The  fragmentation  is  symmetric  across  the  range 
axis  and  asymmetric  across  the  deflection  axis.  The  expected  perforations 
matrices  for  A and  D are  reflections  of  each  other,  and  similarly  for  B 
and  C.  Furthermore,  the  matrices  representing  the  complete  solutions 
(A,C  and  B,D)  are  reflections  of  each  other.  Hence,  the  solutions  for  A 
and  B contain  all  the  information  required  for  a solution  of  either  orien- 
tation over  the  entire  ground  plane.  The  program  has  the  capability  of 
performing  the  symmetry  operations. requi red  to  piece  together  the  matrices 
A and  for  B,  thereby  obtaining  the  complete  solution  for  both  orientations. 
The  output,  including  lethal  area,  average  probability  of  damage,  and 
damage  area,  will  pertain  to  the  complete  ground  plane.  To  use  this 
capability,  we  must  have  0 < IROW  < NXX  (see  paragraph  4.3.2,  card  1 for 
IROW)  for  both  of  the  complete  structure  matrices,  A and  B.  Furthermore, 
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both  values  of  I ROW  should  be  small  because  the  portions  of  the  two 
matrices  A and  B in  rows  numbered  less  than  or  equal  to  the  corresponding 
I ROW  value  are  lost  c the  two  matrices  are  pieced  together.  For  a similar 
reason,  ICOL  should  te  approximately  the  same  for  both  matrices. 

Note  that  this  treatment  depends  upon  the  reflection  symmetry  of 
the  structure  across  the  deflection  axis.  If  the  structure  suffers 
asymmetry  across  that  axis,  a very  similar  treatment  applies,  but  the 
result  has  a different  interpretation.  In  this  event,  structures  B,  D are 
not  actual  structures  with  orientation  opposite  to  A,C.  Rather,  they  are 
hypothetical  structures  which  are  reflections  of  A,C  across  the  range  axis. 
Furthermore,  after  the  matrices  for  A and  B are  pieced  together  at  the 
range  axis,  the  resulting  matrix  A,B  represents  only  A and  C.  If  orientation 
B and  D is  desired,  then  the  solution  for  a second  set  of  conditions 
(structure  B,  reflected  structure  A)  must  be  obtained. 

5.4  Input  Data  Deck. 

This  deck  inputs  a batch  of  cases.  A case  is  one  target  (usually 
consisting  of  more  than  one  component),  projectile,  burst  point  combination. 

A batch  may  contain  as  many  as  25  cases  and,  although  there  is  usually  some 
relationship  among  the  cases  within  a batch,  there  need  not  be.  If  the 
capability  of  paragraph  5.3  is  used,  the  cases  must  come  in  pairs,  and  there 
can  be  25  such  pairs.  In  this  situation,  the  comments  previously  made 
pertaining  to  a case  pertain  to  a case  pair. 

Unless  otherwise  stated,  each  data  field  is  10  columns  wide  and 
each  numeric  datum  is  in  a real  format  with  an  explicit  decimal  point. 

Data  fields  not  specifically  addressed  are  ignored  and  may  be  used  for  any 
desired  purpose. 

5.4.1  The  Data  Deck.  The  data  deck  for  program  Section  II  is  at 
Figure  5-7. 

5.4.2  Definition  of  Each  Input  Variable. 


Card  1: 

Columns  1-10  - the  word  ASSEMBLE. 

KSYM  - KSYM  = 2.  causes  the  union  of  2 expected  perforations 

matrices  at  their  range  axes  to  represent,  over  the  entire 
ground  plane,  a complete  target;  otherwise  ignored. 

NPNTS  - has  two  somewhat  different  interpretations;  if  KSYM  f 2., 

NPNTS  is  the  number  of  cases  in  this  batch,  a case  being  one 
complete  target  at  one  set  of  projectile  conditions;  if 
KSYM  = 2.,  NPNTS  is  the  number  of  case  pairs  in  the  batch; 
in  either  case,  1.  < NPNTS  < 25. 


— 


Columns  21-30  - the  word  CONTOUR;  causes  a contour  map  of  the  probability  of 
damage  matrix  to  be  output  for  each  case  and  damage  criterion 
of  the  current  batch;  also  causes  generation  of  a damage 
function,  storage  thereof  on  the  damage  function  file,  and 
printed  output  thereof;  otherwise,  ignored. 

NSTAPE  - the  unique  whole  number  identifying  the  first  damage  function 

to  be  generated  in  this  batch  of  cases.  The  damage  function 
file  cannot  contain  more  than  5000  records  (functions). 

NSTAPE  can  be  any  whole  number  from  1.  thru  5000.  NSTAPE  is 
the  first  value  of  a consecutive  series  of  numeric  indices 
for  the  functions  to  be  generated  and  filed  in  the  current 
batch  of  cases.  NSTAPE  itself,  however,  need  not  be  consec- 
utive or  even  sequential  with  respect  to  the  functions  stored 
in  the  damage  function  file  by  all  preceding  batches.  If 
NSTAPE  is  not  greater  than  the  index  of  the  highest  indexed 
function  placed  in  the  file  by  all  preceding  batches,  some 
or  all  of  the  earlier  functions  may  be  overlaid  by  functions 
of  the  current  batch.  There  is  a default  value  for  NSTAPE. 

If  NSTAPE  is  blank  or  0.,  its  value  is  computed  to  be  one 
greater  than  the  index  of  the  highest  indexed  function 
already  in  the  file.  If  the  file  is  empty,  the  default  for 
NSTAPE  is  1. 

Columns  51-60  - the  words  NO  PR0BABI ; causes  suppression  of  the  printing  of 
the  probability  of  damage  matrices  for  this  batch;  otherwise, 
ignored. 


Card  2: 

This  is  an  optional  card. 

Columns  1-10  - the  words  EXPECTED  P;  causes  output  of  the  expected  perfor- 

ations matrices  for  individual  cases  of  this  batch  as 
indicated  below: 

Columns  23-30  - in  columns  23-26,  ALL=  and  in  columns  27-30,  a whole  number 
> 0.,  which  number  is  the  power  to  which  10.  will  be  raised 
and  then  used  to  multiply  each  element  of  the  expected  per- 
forations matrix;  the  fractional  part  of  each  element  will 
be  truncated  and  the  matrix  will  be  output  in  16  format. 


Card  3: 

This  card  is  present  only  if  card  2 is  present  and  columns  23-26  thereof 
are  not  as  described  above.  There  must  be  NPNTS  values  of  SH0W(I),  where 
SHOW ( I ) is  defined  in  paragraph  4.3.b.,  card  5;  any  Ith  matrix  for  which 
SH0W(I)  = 0.  will  not  be  output. 


Card  4: 
NKRIT 

KRIT(I) 
Card  5: 


- the  number  of  damage  criteria  which  will  be  considered  for 
each  case  in  this  batch;  L < NKRIT  <_  7. 

- the  Ith  damage  criterion;  KRIT(I)  > 0. 


Cards  5,  6,  and  7 form  a set  which  completely  defines  one  case  in  this 
batch  of  cases.  Either  or  both  of  cards  6 and  7 may  be  omitted  for  any 
card  5,  in  which  case  the  values  from  the  immediately  preceding  card  6 or  7 
are  used.  If  cards  6 or  7 or  both  are  missing  for  the  first  case  of  the 
batch,  the  values  are  taken  to  be  zero.  If  KSYM  = 2.,  a special  rule 
applies.  There  must  be  2 sets  of  cards  5,6,7  to  define  a case;  the  first 
triplet  defines  the  first  half  uf  the  case  pair  and  the  second  triplet,  the 
second  half.  Both  cards  5 must  be  present  for  the  case,  even  if  they  are 
identical  to  each  other.  Cards  6,  7 may  be  present  or  absent  with  the  same 
effect  between  the  halves  of  a case  and  between  cases  as  occurred  between 
cases  for  KSYM  i 2. 


NFILE 


FILE(I) 


Card  6: 
Columns  1-10 
SHFTX(I) 


- the  number  of  logical  records  (expected  perforations  matrices 
for  individual  plates  or  plate  pairs)  required  to  obtain 

the  expected  perforations  matrix  for  the  complete  target 
specified  by  this  set  of  cards  5,  6,  and  7;  1.  £ NFILE  <_  7. 

- the  unique  logical  record  number  associated  with  the 
expected  perforations .matrix  stored  on  the  expected  perfor- 
ations file  for  the  Ith  plate 


plate  (or  plate  pair)  of  the 


current  target. 


- the  words  SHIFT  X 


th 


Card  7: 
Columns  1-10 
SHFTY(I) 


the  number  of  cells  by  which  the  central  cell  of  the  I 
component  plate  (or  plate  pair)  of  this  complete  target  is 
offset  from  the  central  cell  of  the  first  component  of  this 
complete  target  along  the  range  axis;  I must  run  from  2 
to  NPNTS. 


the  words  SHIFT  Y = 

same  as  SHFTX(I),  but  along  the  deflection  axis;  recall 
that  deflection  left  is  positive  and  right,  negative. 


6.  PROGRAM  SECTION  III  - EFFECTIVENESS  SUBMODEL 


6.1  The  DAMAGE  FUNCTION. 

The  damage  functions  generated  by  program  sections  I and  II  are 
not  suitable  for  use  by  this  program  section  because  they  are  not  normal- 
ized, because  their  various  regions  may  not  be  nested  in  a manner  which 
can  be  handled  by  this  program  section,  and  because  they  may  represent 
only  one  half  of  a damage  function.  They  are,  however,  a very  convenient 
point  of  departure  for  the  construction  of  usable  damage  functions. 

There  are  three  types  of  geometry  which  are  allowed  for  the 
regions  of  the  damage  functions.  They  are  identified  as  fully  nested, 
unnested,  and  partially  nested,  and  are  exemplified  in  Figures  6-1,  6-2, 
and  6-3,  respectively.  With  respect  to  Figure  6.3,  note  that  the  damage 
function  must  be  considered  as  having  seven  regions  rather  than  five. 

The  .9  and  yl.  regions  each  must  be  divided  into  two  regions  at  the  range 
axis  so  that  they  can  be  represented  by  rectangles.  There  are  no  restric- 
tions concerning  how  the  probability  of  damage  may  vary  from  one  region 
to  the  next;  however,  there  is  a rule  concerning  geometrical  arrangement 
of  regions.  The  rule  is  that  any  rectangular  contour  which  causes  nesting 
to  occur  must  fully  enclose  all  preceding  contours,  and  each  subsequent 
region  of  that  damage  function  must  completely  enclose  all  regions  pre- 
ceding it.  Once  the  geometry  of  the  final  damage  function  has  been 
determined,  it  must  be  normalized.  That  is,  the  lethal  area  and  damage 
area  for  any  region  of  the  rectangularized,  normalized  damage  function 
must  be  the  same  as  the  sums  of  the  lethal  and  damage  areas,  respectively, 
over  the  elements  of  the  probability  of  damage  matrix  which  were  combined 
to  obtain  that  region.  It  should  be  clear  that  a probability  of  damage 
matrix  or,  more  conveniently,  an  unnormalized  damage  function  from  program 
section  I or  II,  can  result  in  more  than  one  damage  function  as  used  by 
section  III.  Furthermore,  the  choice  of  the  damage  function's  final  form 
is  a matter  of  judgement. 

To  illustrate  the  generation  of  fully  nested  damage  functions, 
we  first  obtain  from  the  unnormalized  damage  function  of  Figure  4-8  what 
is  commonly  called  a "cookie  cutter"  damage  function.  Such  a damage  func- 
tion consists  of  one  region  with  some  constant  probability  of  damage,  and 
the  entire  exterior  region  has  zero  probability  of  damage.  The  total 
lethal  and  damage  areas  are  18.588  and  28.738  square  meters,  respectively, 
so  the  average  probability  of  damage  (abbreviated  PD)  is  .65.  The  only 
quantities  remaining  to  be  defined  are  the  deflection  to  range  ratio  of 
the  rectangle  and  the  location  of  the  rectangle's  center  with  respect  to 
the  burst  location.  We  choose  the  deflection  to  range  ratio  to  be  the 
same  as  it  is  for  the  PD  = .354  region  (3.048/6.096  = .500)  and  we  choose 
the  range  coordinate  of  the  cookie  cutter's  center  to  be  the  same  as  ^or 
the  ?D  = .971  rectangle  (.711  meters  closer  to  the  cannon  than  is  the  burst 
point).  Because  the  target  plate  was  horizontal  and  rectangular,  the 
probability  of  damage  matrix  (Figure  4-7  shows  the  contour  map)  corresponds 
to  a target  enjoying  reflection  symmetry  across  the  range  axis,  and  so 
the  lethal  area,  damage  area,  and  deflection  to  range  ratio  must  be  doubled 
and  the  deflection  coordinate  of  the  cookie  cutter's  center  must  be  set  to 
zero. 
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DEFLECTION 


Figure  6*1.  (U)  One  Half  of  a Fully  Nested  Normalized 
Damage  Function  Enjoying  Reflection 
Symmetry  Across  the  Range  Axis. 


DEFLECTION 


Figure  6-2.  (U)  One  Half  of  an  Unnested  Normalized 

Damage  Function  Enjoying  Reflection 
Symmetry  Across  the  Range  Axis. 


.1 


Figure  6-3.  (U)  A Complete  Partially  Nested  Normalized 
Damage  Function,  Assymetric  Across  the 
Range  Axis. 


RANGE 
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A damage  function  more  closely  approximating  the  probability  of 
damage  matrix  might  be  desired.  In  this  case,  we  might  decide  that  the 
damage  function  should  have  three  regions  corresponding  to  the  PD  = 1.  region, 
the  PD  = .971  through  .449  regions,  and  the  PD  = .354  through  .041  regions. 

We  recognize  at  the  outset  that,  as  above,  the  lethal  area,  damage  area,  and 
deflection  to  range  ratio  for  each  region  will  have  to  be  doubled  and  that 
the  deflection  coordinate  of  the  center  of  each  region  will  be  zero.  The 
innermost  region  (PD  = 1.)  is  easily  done  and  -we  obtain  results  as  shown  in 
Table  6-1, where  all  dimensions  and  coordinates  are  in  meters  and  the  areas 
are  in  square  meters.  The  middle  region  is  done  in  essentially  the  same 
manner,  but  now  six  regions  of  the  unnormalized  damage  function  are  combined. 
We  decide  that  the  deflection  to  range  ratio  should  be  from  the  PD  =..551 
region.  In  obtaining  the  range  and  deflection  dimensions,  we  must  augment 
the  damage  area  by  that  of  the  innermost  region  because  the  three  regions 
are  to  be  completely  nested.  We  choose  the  center  of  the  middle  region  to  be 
at  the  center  of  the  unnormalized  region  with  PD  = .848.  The  outermost 
region  is  obtained  by  combining  the  last  unnormalized  regions,  using  the 
PD  = .146  region  for  both  the  deflection  to  range  ratio  and  the  location  of 
the  region  center. 

The  program  has  the  capability  of  performing  the  preceding  com- 
putations (paragraph  6.4.2,  card  14).  The  resulting  normalized  damage 
function  will  be  fully  nested.  The  input  (unnormalized)  damage  function  is 
read  in  from  the  damage  function  file  and  the  regions  thereon  are  assumed 
to  be  ordered  from  outermost  to  innermost,  which  is  the  sequence  in  which 
they  were  written  on  the  damage  function  file  by  program  sections  I and  II. 

To  obtain  the  three  region  damage  function,  we  use  STEPS 2(1)  = 1.,  STEPS2 ( 2 ) 

= 6.,  and  STEPS2(3)  = 4.  Note  that  the  elements  of  the  STEPS2  vector  are 
ordered  from  innermost  to  outermost  (usually  high  PD  to  low  PD)  regions  of 
the  desired  normalized  damage  function  and  that  their  sum  is  equal  to  the 
number  of  regions  in  the  unnormalized  damage  function  on  file.  If  the  sum 
over  the  elements  of  the  STEPS2  vector  is  less  than  the  number  of  regions 
in  the  function  on  file,  then  the  innermost  N = STEPS2(J)  reaions  in  the 

function  on  file  will  be  used  and  the  remainder  ignored. 

In  either  case,  the  program  will  combine  regions  as  indicated  by 
the  STEPS2  vector  and  will  obtain  the  deflection  to  range  ratio  for  each 
resulting  region  from  the  outer  boundaries  of  that  region.  Normalization 
will  occur.  The  midpoint  of  each  region  will  be  computed  as  the  midpoint 
of  each  of  the  normalized  rectangles.  If  the  "symmetry"  option  (paragraph 
6.4.2,  card  14)  is  used,  the  lethal  area,  damage  area,  and  deflection  to 
range  ratio  will  be  doubled  and  the  deflection  coordinate  of  each  region's 
midpoint  will  be  set  equal  to  zero.  Working  from  innermost  to  outermost, 
the  program  will  compare  the  (i  + l)™  region  to  the  i*-'1  region  to  insure 
that  all  regions  are  nested.  If  the  (i  + l)th  rectangle  cannot  enclose 
the  ith  rectangle  because  of  shape,  the  shape  of  the  (i  + l)t!l  rectangle  will 
be  altered  just  enough,  maintaining  normalization,  to  allow  it  to  enclose 
the  ith  rectanale.  Also,  if  the  (i  + l)th  rectangle  must  be  shifted  to 
enclose  the  itff  rectangle,  the  (i  + l)th  midpoint  will  be  shifted  just 
enough  for  nesting  to  occur. 


TABLE  6-1.  A FULLY  NESTED  DAMAGE  FUNCTION 
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For  unnested  damage  functions,  and  for  the  unnested  part  of  a 
partially  nested  damage  function,  the  procedure  is  similar,  but  less 
complicated  simply  because  we  need  not  be  concerned  about  nesting.  How- 
ever, the  program  has  no  capability  to  generate  such  damage  functions. 

In  any  event,  the  program  reflects  the  complete  damage  function  in 
its  final  form  through  the  projectile  burst  point.  This  transforms  the 
function  from  one  describing  the  damage  sustained  by  a target  for  every 
location  of  that  target  about  the  burst  point  into  a function  describing 
the  damage  sustained  by  a target  due  to  a projectile  burst  at  any  point 
about  that  target,  with  constant  HOB. 

6.2  The  Program. 

This  section  of  the  program  uses  the  methodology  of  reference  4, 
with  some  changes  in  the  computer  program  to  increase  flexibility  and  delete 
unneeded  portions,  to  compute  the  probability  that  each  target  in  a group  of 
several  similar  targets  is  damaged.  This  probability  is  stored  internally 
as  an  element  of  a 3-dimensional  matrix,  henceforth  called  the  S matrix. 

The  three  subscripts  of  these  matrix  elements  identify  the  target  element 
sustaining  damage,  the  size  of  the  volley  which  caused  the  damage,  and  the 
HOB  of  all  projectiles  in  the  volley. 

Essentially,  the  S matrix  is  computed  as  follows.  The  specified 
number  of  projectiles  is  fired  at  the  selected  aim  points;  proper  launch 
is  assumed.  Round- to-round  and  occasion-to-occasion  er-or  (assumed  normal) 
and  projectile  functioning  reliability  (binomial)  distributions  are  applied 
to  the  projectiles.  HOB  is  constant.  A distribution  of  burst  points  over 
the  specified  target  array  is  obtained.  The  burst  point  distribution  and 
the  damage  function  representing  the  target  element  are  combined  and 
integrated  over  the  ground  plane.  The  result  is  one  element  of  the  S matrix 
Results  are  obtained  for  each  volley  size,  each  target  in  the  target  array, 
and  each  HOB  (each  HOB  has  its  own  damage  function)  to  obtain  all  elements 
of  the  S matrix. 

Elements  of  the  S matrix  can  be  averaged  over  any  specified  HOB 
distribution  (so  long  as  the  distribution  includes  only  HOB's  for  which 
there  are  elements)  to  obtain  a two-dimensional  S'  matrix.  The  two  sub- 
scripts of  these  matrix  elements  identify  the  target  element  sustaining 
damage  and  the  size  of  the  volley  causing  the  damage. 

Optimization  can  be  performed  upon  the  S matrix  for  a specified 
volley  size  to  determine  which  HOB  causes  the  most  damage.  This  HOB  is 
then  treated  as  an  HOB  distribution  with  standard  deviation  equal  to  zero 
to  obtain  an  S'  matrix.  The  optimization  consists  of  summing  the  S matrix 
elements  over  the  various  targets  for  each  HOB  at  the  specified  volley  size. 
The  HOB  yielding  the  largest  sum  is  defined  as  optimum  for  that  volley  size. 

From  the  S'  matrix  are  obtained,  for  each  volley  size  specified, 
the  probabilities  that  various  numbers  of  targets  in  the  target  array  sus- 
tain damage,  the  average  probability  (averaged  over  the  target  elements) 


54 


that  some  target  element  sustains  damage,  and  the  most  and  least  damageable 
target  elements  with  the  probability  that  they  sustain  damage. 

Consider  the  plane  of  elements  in  the  S matrix  corresponding  to 
one  target  element.  Each  such  plane  is  generated  by  a machine  computation 
which  can  be  time  consuming.  Under  certain  conditions,  the  number  of  such 
planes  to  be  computed  can  be  halved  and  the  remainder  obtained  through  an 
appeal  to  reflection  symmetry.  The  conditions  are  that  the  damage  functions 
for  all  specified  HOB's,  the  aim  point  array,  and  the  target  array  enjoy 
reflection  symmetry  across  the  range  axis.  Note  that,  apart  from  being 
common  to  both  the  aim  points  and  the  target  elements,  the  Cartesian 
coordinate  system  used  is  arbitrary  and  its  origin  may  be  chosen  in  any 
convenient  manner. 

When  these  conditions  are  satisfied,  the  complete  aim  point  array 
is  input,  but  only  one  of  the  two  symmetric  halves  of  the  target  array  is 
input,  and  the  SYMMETRIC  ARRAY  option  of  input  data  cards  16  and  19  through 
21  (Figure  6-4)  is  used.  This  allows  the  use  of  a target  array  containing 
twice  as  many  target  elements  (50)  as  can  be  input  (25),  and  results  in 
some  economy  of  machine  time. 

6.3  Output. 

The  input  data  are  printed  out.  The  S and  S'  matrices  are  printed 
out.  Depending  upon  the  number  of  targets  in  the  target  array,  such  of  the 
following  as  apply  are  printed  out  for  all  specified  volley  sizes  at  the 
specified  HOB  distribution  (the  optimum  HOB,  if  computed,  being  considered 
as  an  HOB  distribution  with  zero  standard  deviation).  These  and  the  S' 
matrix  are  printed  out  only  when  requested  by  the  appropriate  cards  16 
through  21  (see  paragraph  6.4.2),  whereas  the  S matrix  is  always  printed  out. 

a.  Probability  that  no  targets  are  damaged. 

b.  Probability  that  exactly  one  target  is  damaged. 

c.  Probability  that  exactly  two  targets  are  damaged. 

d.  Probability  that  at  least  one  target  is  damaged. 

e.  Probability  that  at  least  two  targets  are  damaged. 

f.  Probability  that  at  least  three  targets  are  damaged. 

g.  Probability  that  less  than  two  targets  are  damaged. 

h.  Probability  that  less  than  three  targets  are  damaged. 

i.  Average  probability  that  a target  is  damaged. 

j.  Probability  that  the  most  vulnerable  target  is  damaged. 

k.  Probability  that  the  least  vulnerable  target  is  damaged. 


6 . 4 Input  Data  Deck. 

This  deck  inputs  a batch  of  cases.  Data  cards  1 through  15  provide 
data  used  by  all  cases  in  the  batch.  These  data  include  the  number  of  cannon 
and  their  aim  points,  the  number  of  targets  in  the  target  array  and  their 
locations,  the  number  of  rounds  fired  by  each  cannon  for  each  of  the 
volley  sizes,  ballistic  data,  functioning  reliability,  and  damage  functions 
for  each  of  as  many  as  25  HOB's.  Data  cards  16  through  18  or  cards  19  or  20 
or  21  define  a specific  case  within  the  batch.  As  many  cases  as  desired 
may  be  run  within  one  batch.  A case  now  means  a specific  HOB  distribution 
consisting  of  some  or  all  of  the  HOB's  for  which  damage  functions  were  input. 

Unless  otherwise  stated,  each  data  field  is  10  columns  wide,  and 
each  numeric  datum  is  in  a real  format  with  an  explicit  decimal  point.  Data 
fields  not  specifically  addressed  are  ignored  and  may  be  used  for  any  desired 
purpose. 


Figure  6-4. 


6.4.1  The  Data  Deck.  The  data  deck  for  program  section  III  is  at 


6.4.2  Definition  of  Each  Input  Variable. 


Card  1: 


Columns  1-10  - the  word  ENCOUNTER 

N - the  number  of  cannon  to  fire;  1.  < N < 100. 

NVEH  - the  number  of  targets  in  the  target  array  to  be  fired  upon; 

1.  <_  NVEH  ^25.;  see  card  16,  columns  31-40. 

NVS  - the  number  of  volley  sizes  to  be  solved;  1.  £ NVS  <_  25. 

NHTS  - the  number  of  HOB's,  and  hence,  the  number  of  damage  functions 

to  be  solved;  1.  £ NHTS  <_  25. 

U1  - the  number  of  increments  for  the  range  and  deflection  integrals 

over  a quarter  of  the  integration  plane;  usually,  U1  > 8. 

Cards  2 and  3: 

A(I),  B ( I ) - the  range  and  deflection  coordinates,  respectively,  of  the 
aiming  point  for  the  Ith  cannon;  in  meters;  recall  that 
deflection  left  is  positive  and  right,  negative. 

There  may  be  as  many  as  13  cards  for  the  A's  and  the  same  number  for  the  B's 
to  allow  the  input  of  N values  of  A and  B. 

Card  4; 

An  optional  card;  has  the  words  SAME  CANN0  in  columns  1-10.  If  present, 
cards  2 and  3 must  be  absent.  Causes  the  aiming  points  of  the  immediately 
preceding  section  III  data  deck  to  be  used. 


FIGURE  6-<t.  DATA  DECK  FOR  PROGRAM  SECTION  III 


Cards  5 and  6: 

VOSR(I),  VOSD(I)  - the  range  and  deflection  coordinates,  respecti velv , of 
the  I***1  target  in  the  target  array;  in  meters;  recall 
that  deflection  left  is  positive  and  right,  negative. 

There  may  be  as  many  as  four  cards  for  the  VOSR's  and  the  same  number  for 
the  VOSD's  to  allow  for  the  input  of  NVEH  values  of  VOSR  and  VOSD.  See 
card  16,  columns  31-40. 

Card  7: 

An  optional  card;  has  the  words  SAME  TARGE  in  columns  1-10.  If  present, 
cards  5 and  6 must  be  absent.  Causes  the  target  locations  of  the  immediately 
preceding  section  III  data  deck  to  be  used. 

Card  8: 

NN(I)  - the  number  of  rounds  to  be  fired  by  each  cannon  in  the 

I^  volley  size  to  be  solved. 

There  may  be  as  many  as  four  cards  for  the  NN's  to  allow  for  the  input  of 
NVS  volley  sizes. 

Card  9: 

An  optional  card;  has  the  words  SAME  VOLLE  in  columns  1-10;  if  present, 
card  8 must  be  absent.  Causes  the  volley  sizes  of  the  immediately  preceding 
section  III  data  deck  to  be  used. 

Card  10: 

SI,  S2 


Tl,  T2 


R 


Card  11: 

An  optional  card.  If  present,  card  10  must  be  absent;  has  the  words  SAME 
| ERROR  in  columns  1-10.  Causes  SI,  S2,  Tl,  T2  and  R from  the  immediately 

preceding  section  III  data  deck  to  be  used. 
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- precision  errors;  the  range  and  deflection  standard 
deviations  of  the  distribution  of  one  cannon's  burst 
points  at  one  HOB  on  one  occasion  about  their  own  mean; 
in  meters. 

- mean  point  of  impact  (MPI)  errors;  the  range  and  deflec- 
tion standard  deviations  of  the  distribution  of  one 
cannon's  single  occasion  mean  burst  points  at  one  HOB 
about  their  own  mean;  in  meters. 

- the  complete  round  functioning  reliability,  given  proper 
launch. 


Card  12: 

Cards  12  and  13  define  one  damage  function  in  final  form.  There  must  be 
NHTS  sets  of  cards  12  and  13.  These  sets  need  not  be  in  any  particular 
order. 


HT(  I ) 


NUMB ( I ) 


NUMNN(I) 


Card  13: 
XSSP(I,J) 


- the  HOB  to  which  the  I damage  function  applies;  any 
two  elements  must  differ  from  each  other  by  at  least  .1% 
of  their  sum;  in  meters. 

- the  number  of  regions  in  the  Ith  damage  function;  1. <_ 
NUMB(I)  < 13. 

th 

- the  number  of  unnested  regions  in  the  I damage  function; 
0.  < NUMNN(I)  < NUMB ( I ) . 


XSSP(I,J)  - the  average  probability  of  damage  over  the  JUI  region  of 

the  Itf1  damage  function. 

XBB3( I ,J) , XBB4 ( I ,J ) - the  range  and  deflection  dimensions , respecti vely , 
of  the  region  of  the  I tn  damage  function;  in  meters. 

XR0SE(I,J),  XDOSE ( I , J ) - actually  XROSET(I.J),  XDOSET(I.J);  the  range  and 
deflection  offsets,  respectively,  of  the  center  of  the 
region  of  the  I™  damage  function  with  respect  to  the 
target  location;  in  meters. 

There  must  be  NUMB(I)  cards  13  in  the  1^  set  of  cards  12  and  13.  Their 
order  must  be  from  innermost  to  outermost  damage  function  region.  The  order 
among  themselves  of  those  cards  defining  unnested  regions  is  unimportant, 
but  all  such  unnested  cards  must  occur  before  any  cards  defining  nested 
regions  in  the  lzn  set. 

Card  14: 

One  card  14  replaces  one  set  of  cards  12  and  13.  Any  and  any  number  of 
those  card  sets  may  be  so  replaced.  Causes  input  of  data  from  the  damage 
function  file,  from  which  data  the  program  will  generate  a fully  nested 
damage  function. 


the  Jtn  region  of 


Columns  1-10 


Columns  11-20 


any  alphanumeric  character  string  which  cannot  be 
interpreted  as  a number. 

if  columns  11-20  contain  the  word  SYMMETRY,  the  generated 
damage  function  will  be  assumed  to  be  half  of  a damage 
function  enjoying  reflection  symmetry  across  the  range 
axis  and  it  will  be  reflected  across  the  range  axis  to 
produce  a complete  function.  Lethal  area  and  damage  area 
will  double.  Otherwise,  ignored. 
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FILE2 


ST  (J ) 


5 column  whole  number;  the  unique  logical  record  number 
identifying  the  desired  d.-'ta  on  the  damage  function  file. 

5 column  whole  number;  actual ly  STEPS2(J);  one  value  for 
each  region  desired  in  the  normalized  fully  nested  damage 
function  to  be  generated;  each  value  indicates  the 
adjacent  regions  of  the  unnormalized  damage  function  on 
file  which  will  be  combined  to  arrive  at  the  Jth  region 
of  the  desired  normalized  damage  function;  the  elements 
of  this  set  must  be  ordered  from  innermost  to  outermost 
regions  of  the  desired  function;  the  sum  of  all  the 
STEPS2(J)  must  not  be  greater  than  the  number  of  regions 
in  the  function  on  file  (not  more  than  11  regions);  if 
the  sum  of  all  the  STEPS2(J)  is  less  than  the  number  of 
regions  in  the  function  on  file,  then  only  the  innermost 
N = j STEPS2(J)  regions  of  the  damage  function  on  file 
will  be  used;  the  first  STEPS2 ( J ) = 0.  or  blank  terminates 
the  vector.  There  is  one  exception  to  the  preceding 
rules;  if  aJJ_  the  STEPS2(J)  are  0.  (or  blank)  then  no 
regions  on  file  will  be  combined;  that  is,  the  result  is 
identical  with  the  STEPS2  vector  being  all  l.'s  and  having 
as  many  elements  as  there  are  regions  in  the  damage 
function  on  file. 


Card  15: 

An  optional  card.  If  present,  all  cards  12,  13,  and  14  must  be  absent. 

The  words  SAME  DAMAGE  in  columns  1-10  cause  the  entire  set  of  normalized 
damage  functions  from  the  immediately  preceding  section  III  data  deck  to  be 
used. 

Card  16: 


Cards  16,  17,  18  form  a set  defining  one  HOB  distribution.  There  may  be  as 
many  or  as  few  such  sets  as  desired,  or  none. 


Columns  1-10 


Columns  31-40 


- the  words  HOB  DISTRI. 

- the  number  of  HOB's  in  the  distribution;  1.  < NH0B  < 25. 

- the  word  SYMMETRIC;  causes  the  target  array  to  be  reflec- 
ted across  the  range  axis,  thus  producing  a new  target 
array  enjoying  reflection  symmetry  across  the  range 

axis  and  having  twice  as  many  target  elements  as  the 
array  which  was  input;  can  be  used  only  when  the  aim- 
points  and  the  damage  function  are  symmetric  across  the 
range  axis;  does  not  carry  over  from  one  case  to  the  next 
or  from  one  batch  to  the  next;  otherwise,  ignored. 
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Columns  51-60  - the  word  OPTIMIZE;  causes  HOB  optimization  without  prior 
use  of  an  HOB  distribution. 


L - the  number  of  rounds  fired  per  cannon  at  which  optimization  will 

be  performed;  if  L is  not  a member  of  the  NN ( I ) vector,  NN ( 1 ) will 
be  used;  ignored  if  columns  51-60  are  ignored. 

Cards  17  and  18: 

HOB (I ) , DF ( I ) - the  Ith  HOB,  and  its  weighting  factor,  respectively,  for 
the  current  HOB  distribution  function;  each  HOB  must  be 
within  .1%  of  an  element  of  the  HT(I)  vector. 

There  may  be  as  many  as  four  cards  17  and  four  cards  18  to  allow  the  input 
of  NHOB  values  of  HOB  and  corresponding  values  of  DF. 

Card  19: 

Cards  19,  20,  and  21  each  perform  a function  equivalent  to  that  of  the  set 
of  cards  16,  17,  and  18.  As  many  or  as  few  as  desired  may  be  used  and  they 
may  be  intermingled  with  each  other  and  with  the  set  of  cards  16,  17,  and  18, 
in  any  order  desired  (but  not  inside  the  set  16,  17,  and  18),  with  the 
obvious  restriction  on  card  20,  that  there  must  have  been  an  immediately 
preceding  case. 

Columns  1-10  - the  words,  ALL  HOBS;  causes  the  current  HOB  distribution  to 
be  a uniform  distribution  consisting  of  all  HOB's  in  the  set 
(HT). 

Columns  31-80  - same  as  card  16. 

Card  20: 

Columns  1-10  - the  words,  SAME  HOB  D;  causes  the  current  HOB  distribution 
to  be  the  same  as  in  the  immediately  preedino  case;  the 
immediately  preceding  case  can  be  the  last  i-.  e of  the 
immediately  preceding  batch. 

Columns  31-80  - same  as  card  16. 

Card  21: 

Columns  1-10  - the  word,  OPTIMIZE;  causes  HOB  optimization  without  prior 
use  of  an  HOB  distribution. 

L - same  as  card  16. 

Columns  31-50  - same  as  card  16. 
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7.  PROGRAM  SECTION  IV  - PERFORATIONS;  POSSIBILITY  AND  SIGNIFICANCE 


7.1  The  Program. 

This  program  section  provides  two  functions  not  directly  related 
to  those  of  program  sections  I,  II,  and  III.  It  predicts,  for  an  armor 
plate  and  a fragment  distribution,  whether  perforations  can  ever  occur  except 
for  direct  hits  and  very  near  misses.  If  perforations  can  occur,  it  computes 
the  average  probability  of  incapacitation  of  a soldier  protected  against  the 
fragment  distribution  by  the  armor  plate,  assuming  that  the  soldier  will  be 
hit  by  the  distribution  of  residual  fragments.  Once  again,  the  methodology 
is  that  of  Thor  Reports  47  and  70  (references  1 and  2).  Air  drag  is 
neglected.  For  those  fragments  which  perforate,  the  residual  mass  and  vel- 
ocity are  computed.  For  each  perforating  fragment,  the  methodology  of 
reference  5 is  used  to  obtain  the  probability  of  incapacitation  of  a soldier, 
given  that  he  is  hit  by  the  residual  fragment.  This  methodology  is  augmented 
by  the  computation  of  a V™  using  the  equations  of  reference  8.  Correspon- 
dence of  the  naked  soldieruand  the  soldier  wearing  a helmet  and  winter 
uniform  in  reference  5 to  the  soldier  wearing  a summer  uniform  and  the 
soldier  wearing  a winter  uniform  in  reference  8,  respectively,  is  assumed. 
These  P.,  's  are  averaged  over  the  distribution  of  perforating  fragments  to 
obtain  /Hthe  average  PT/H  per  fragment  for  the  soldier  protected  by  the 
armor  plate. 

7.2  Output. 

The  input  data  are  printed  out.  Computed  output  is  shown  in 
Figure  7-1. 

7.3  Input  Data  Deck. 

This  data  deck  is  identified  by  its  first  card  which  contains  the 
words  SECTION  IV  in  columns  1-10.  After  this  card  are  input  any  number  of 
"batches"  of  cases.  One  batch  is  terminated  by  a card  with  the  word 
COMPUTE  in  columns  1-10.  One  case  is  one  armor  plate  or  plate  assembly  at 
one  obliquity  angle.  For  any  given  batch,  the  input  data  are  divided  into 
six  data  card  subsets,  which  are  the  fragmentat.on  data,  terminal  velocity, 
Thor  constants,  angles,  armor,  and  incapacitation  constants  subsets.  These 
subsets  may  be  presented  in  any  order. 

All  six  subsets  must  be  present  for  the  first  batch,  except  that, 
if  the  desired  fragmentation  data  are  currently  in  storage  from  program 
section  I (see  paragraph  8),  then  the  fragmentation  data  subset  need  not  be 
present.  For  subsequent  batches,  only  those  data  card  subsets  whose  input 
data  are  to  be  changed  need  be  present. 

Unless  otherwise  stated,  each  data  field  is  10  columns  wide,  and 
each  numeric  datum  is  in  a real  format  with  an  explicit  decimal  point. 

Data  fields  not  specifically  addressed  are  ignored  and  may  be  used  for  any 
desired  purpose. 

7.3.1  The  Data  Deck . The  data  deck  for  program  section  IV  is  at 

Figure  7-2. 
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figure  7-1.  sample  output  gf  program  section  iv 


mm 


7.3.2  Definition  of  Each  Input  Variable. 

Card  1: 

Columns  1-10  - the  words  SECTION  IV. 

The  following  card  sets  define  a batch  of  cases.  The  card  sets  may  appear 
in  any  order.  The  batch  is  terminated  by  a COMPUTE  card.  As  many  batches 
as  desired  may  follow  the  SECTION  IV  card. 

Card  set  A: 

This  card  set  inputs  the  fragmentation  data  and  is  identical  with 
cards  11  through  14  of  paragraph  4.3.2. 

Card  set  B: 

Columns  1-10  - the  words  TERMINAL  V. 

VTERM  - the  projectile  terminal  velocity  in  feet  per  second. 

Card  set  C: 

This  card  set  inputs  as  many  as  20  sets  of  Thor  constants  from 
which  are  drawn  the  appropriate  sets  for  the  various  armor  plates  or  armor 
plate  assemblies  in  this  batch. 

Card  Ci: 

Columns  1-10  - the  words  THOR  DATA. 

NTHOR  - the  number  of  sets  of  Thor  constants;  1.  < NTHOR  < 20. 


Cards  C2  and  C3: 

This  card  pair  defines  one  set  of  Thor  constants.  The  cards  are 
identical  with  cards  8 and  9 of  paragraph  4.3.2.  There  must  be  NTHOR  such 
card  pairs. 

Card  set  D: 

i 

Columns  1-10  - the  word  ANGLES. 

NTHETA  - the  number  of  obliquity  angles  to  be  considered  for  each  armor 
plate  or  plate  assembly;  1.  £ NTHETA  £ 20. 


, - THETA(I)  - the  I obliquity  angle  to  be  considered  for  each  armor  plate 

» or  plate  assembly;  in  degrees;  there  must  be  NTHETA  such 

! anal es : 0.°  < THETA  < 90°. 


This  card  set  defines  the  armor  plates  processed  in  this  batch. 


Card  El: 


Columns  1-10  - the  word  ARMOR. 
NPLATE  - the  number  of  plates. 


Card  E2: 


This  card  defines  an  armor  plate.  Exactly  NPLATE  such  plates  must 
be  defined. 

PLATE(I)  - 10  alphanumeric  characters;  specifies  the  material  of 
which  the  I™  plate  is  made;  must  match,  character  by  character,  one  of  the 
variables  METAL(K)  from  card  set  C. 

T ( I ) - the  thickness,  in  inches,  of  the  plate. 

Card  set  F: 

This  card  set  defines  the  incapacitation  parameters  used  in  computing  Pj/H- 
Columns  1-10  - the  word  INCAPACITA. 

Columns  31-40  - the  wore  SUMMER  or  the  word  WINTER;  describes  the 
type  of  uniform  being  worn. 

C0EF1,  C0EF2,  C0EF3  - three  constants  used  to  compute  PT/H  for 
the  selected  kill  criterion  (see  reference  5 and  appendix  B). 

Last  card  of  any  batch: 

Columns  1-10  - the  word  COMPUTE. 


8.  THE  COMPLETE  JOB  DECK 
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8. 1 The  Data  Deck. 

Thus  far,  the  four  sections  of  the  program  and  their  data  decks 
(herein  called  subdecks)  have  been  discussed  separately.  In  an  actual 
machine  run,  any  number  of  each  of  the  four  types  of  data  subdecks  may  be 
stacked  in  any  order  desired,  so  long  as  one  insures  that,  when  a subdeck 
is  processed,  it  will  not  demand  information  produced  by  some  later  subdeck. 

Each  type  of  subdeck  has  certain  options  allowing  the  use  of  data 
read  in  from  a preceding  subdeck  of  the  same  type.  An  important  constraint 
upon  the  complete  data  deck  is  that  each  of  the  four  program  sections 
destroys  the  input  data  belonging  to  any  previous  different  programTection, 
with  one  exception,  which  is  discussed  below. 

Program  sections  I and  IV  do  not  destroy  each  other's  static  frag- 
mentation data.  That  is,  either  of  these  program  sections  can  use  static 
fragmentation  data  read  in  by  itself  or  by  the  other  and  will  leave  that 
data  ready  to  be  used  again  by  itself  or  by  the  other.  Any  time  either  of 
these  program  sections  reads  in  such  data,  it  replaces  the  fragmentation  data 
previously  read  in  by  either  section. 

A complete  data  deck  must  contain  at  its  end  either  of  two  cards 
not  heretofore  discussed. 

One  of  these  is  a card  with  the  words  END  OF  DAT  in  columns  1 
through  10.  This  card  terminates  the  data  deck.  Succeeding  cards  are 
ignored. 


The  other  card  can  have  any  of  the  forms  shown  in  Figure  8-1. 
Only  columns  1-10  and  41-50  are  read;  other  columns  are  ignored.  These 
cards  cause  termination  of  the  run  immediately  after  completion  of  the 
specified  operations. 

EXPECTED  P causes  an  index  of  the  expected  perforations  matrix 
file  to  be  printed.  An  example  is  at  Figure  C-13.  The  column  headings 
have  the  following  meanings: 

FILES  is  the  unique  integer  identifying  an  expected 
perforations  matrix;  used  as  an  element  of  the  set  FILE 
(paragraph  5.4.2,  card  5). 

NXX,  NYY  are  as  in  paragraph  4.3.2,  card  1. 

CELL  is  the  edge  length  of  the  cell  with  which  the  ground 
plane  is  gridded,  in  inches. 

HOB  is  the  height  of  burst  above  the  ground,  in  meters. 

I ROW,  ICOL  are  as  in  paragraph  4.3.2,  card  1. 
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LYMYB,  LYMYT,  LYMXL,  LYMXR  are  the  row  and  column  numbers 
which  define  the  minimum  and  maximum  deflection  and  range 
limits,  respectively,  of  the  non-zero  portion  of  the 
expected  perforations  matrix.  If  the  matrix  is  all  zeroes, 
then  these  numbers  merely  define  a small  region  about 
the  projectile. 

A,  B,  NSIDE,  NX,  NY,  NZ,  T are  elements  of  the  sets  ALPHA, 
BETA,  NSIDE,  NXTGT,  NYTGT,  NZTGT , and  TT,  respectively 
(paragraph  4.3.2,  card  5). 

THCK  is  an  element  of  the  set  TTHCK  (paragraph  4.3.2,  card 
3). 

DAMAGE  FUN  causes  the  contents  of  the  damage  function  file  to  be 
printed.  An  example  is  at  Figure  C-14.  The  column  headings  have  the 
following  meanings: 

FILE  NUMBER  is  the  unique  integer  identifying  a damage 
function;  used  as  FILE2  (paragraph  6.4.2,  card  14). 

POINT  NUMBER  identifies  the  Ith  case  of  a batch  of  cases 
from  program  section  II. 

SYMMETRY  CHARACTER  is  KSYM  (paragraph  6.4.2,  card  1). 

DAMAGE  CRITERION  is  the  number  of  perforations  defined  as 
the  damage  criterion. 

LYMXL,  LYMXR,  LYMYB,  LYMYT  are  as  above. 

PROJECTILE  LOCATION  is  the  range,  deflection  coordinate 
pair  of  the  burst  point  with  respect  to  the  probability 
of  damage  function. 

HOB  is  the  height  of  burst  above  the  ground,  in  meters. 

CELL  is  as  above. 

NUMBER  OF  REGIONS  is  the  number  of  different  valued 
probability  of  damage  regions  in  the  probability  of  damage 
function . 

THE  DAMAGE  FUNCTION  REGIONS  data  completely  defines  the 
unnormalized  damage  function. 

8.2  The  Control  Deck. 

Certain  control  cards  are  needed  to  cause  problem  loading  and 
execution  and  to  insure  availability  and  permanence  of  files.  Depending 
on  the  situation,  these  functions  can  be  accomplished  in  any  number  of 


'*»**•* 


ways,  and  hence  will  not  be  discussed  in  detail.  The  control  decks  used 
to  produce  the  example  of  Appendix  C are  presented  and  explained  in 
Appendix  D.  One  point,  however,  is  worthy  of  note.  On  several  figures 
providing  examples  of  output  in  the  body  of  the  report  and  in  Appendix  C, 
the  names  assigned  by  the  user  to  his  perforation  matrix  and  damage  func- 
tion files  appear.  These  names  are  made  available  to  the  program  by  the 
control  deck.  Note  from  Appendix  A that  the  first  FORTRAN  statement  in 
the  program  is: 

PROGRAM  LASVEM  (INPUT,  OUTPUT,  TAPE1 , TAPE2 , . . . ) 


where: 


INPUT  is  the  system  input  file. 

OUTPUT  is  the  system  output  file. 

TAPE1  is  the  perforation  matrix  file. 

TAPE2  is  the  damage  function  file. 

The  user's  local  file  names  for  the  perforation  matrix  file  and  the  damage 
function  file  must  be  equi valenced, by  appropriate  commands,  to  TAPE1  and 
TAPE2,  respectively. 
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APPENDIX  B 

PROBABILITY  OF  INCAPACITATION  EQUATIONS 

This  appendix  presents  sufficient  information  concerning  the 
probability  of  incapacitation  equations  for  soldiers  (references  5 and  8) 
to  allow  those  equations  to  be  encoded  in  a manner  compatible  with  the 
LASVEM  computer  program.  The  equations  themselves  have  been  deleted  from 
this  report  in  order  to  keep  it  unclassified. 

Cards  PKFCN  6 with  cards  PKFCN  20,  PKFCN  21,  PKFCN  23,  and 
PKFCN  24  compute  a V5Q  level  as  described  in  paragraph  1 of  reference  8. 
That  reference/  uses  a factor,  'A',  the  mean  presented  area  in  centi- 
meters squared.  Centimeters  must  be  converted  to  inches  and  the 
presented  area  must  be  eliminated: 


■(: 


M_ 

0610234 


SHAPE, 


Equation  B-l 


Note  that  card  SEXl|/l29  incorporates  the  factor,  .0610234,  and  that  card 
SEXU188,  which  is  the  call  to  function  PKFCN,  defines  AA,  BB  and,  in  turn 
A,  B,  as  fragment  mass  and  speed,  respectively. 

Cards  PKFCN  8 and  PKFCN  9 compute  the  probability  of  incapacita- 
ting a soldier  with  a fragment  of  mass  A (in  grains)  and  velocity  B (in 
fps).  Reference  6 uses  the  three  variables  a,  b,  and  n which  are 
encoded  as  C0EF1 , C0EF2,  and  C0EF3,  respectively. 
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APPENDIX  C 


AN  ILLUSTRATIVE  EXAMPLE 


This  appendix  provides  an  example  of  the  use  of  this 
methodology  from  start  to  finish.  Although  the  example  is  typical, 
it  is  very  simple.  None  of  the  data  necessarily  represent  actual 
or  developmental  systems.  Not  all  the  options  of  the  program  are 
exercised,  nor  is  all  of  the  output  shown. 

We  wish  to  determine  measures  of  the  damage  sustained  by  a 
force  of  four  APC's  when  the  force  is  fired  upon  by  artillery  firing 
HE  shell  fuzed  for  air  burst.  The  APC's  are  oriented  in  the  direction 
of  fire.  A scenario  such  as  this  might  be  useful  for  determining  the 
protection  provided  to  the  occupants  by  the  APC  armor  against  friendly 
supporting  fires.  The  scenario  is  depicted  in  Figure  C-l. 

Various  numbers  of  rounds  will  be  fired  by  one  howitzer 
battery  consisting  of  six  cannon.  The  aiming  policy  is  battery  parallel 
sheaf  with  a hypothetical  cannon  at  the  battery's  geometrical  center 
aiming  at  the  center  of  the  APC  force.  Ballistic  data  and  fuze  data 
are  shown  in  Table  C-l.  APC  structure  is  indicated  by  a shop  diagram 
at  Figure  C-2.  Data  for  the  individual  APC  armor  plates  are  in 
Table  C-2.  We  use  one  perforation  as  the  damage  criterion. 

First,  program  section  IV  will  be  used  to  determine  which 
plates  have  any  chance  of  sustaining  perforations.  There  are  two 
reasons  for  using  section  IV: 

1.  Plates  identified  as  not  sustaining  any  perforations  can 
be  eliminated  from  consideration,  thus  conserving  both  machine  time  and 
user  effort. 

2.  An  estimate  of  the  average  probability  of  incapacitation 
of  soldiers  hit  by  residual  fragments  will  be  produced. 

The  input  data  deck  is  at  Figures  C-3  and  C-4.  Plates  D and 
E are  not  considered  at  all  because  they  are  shielded  from  fragmentation 
by  suspension  and  mobility  components  and  from  direct  hits  by  the  rest 
of  the  APC.  Data  from  other  sources  indicates  that  the  spaced  armor  will 
not  sustain  perforations.  Results  of  this  first  step  are  summarized  in 
Table  C-3.  A sample  of  the  output  is  at  Figure  C-5. 

Plates  C,  F,  and  H sustain  no  damage  from  fragmentation; 
however,  they  may  sustain  massive  damage  from  direct  hits  and  very  near 
misses  (see  paragraph  4.1).  Although  we  may  process  all  three  of  these 
plates  for  direct  hits  and  very  near  misses  without  error,  we  need 
process  only  plate  F,  because  the  hits  and  very  near  misses  on  plate  C 
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TABLE  C-l.  BALLISTIC  AND  FUZE  DATA 


Projectile  Terminal  Velocity 

Projectile  Angle  of  Fall 

Fuze  Reliability 

Fuze  HOB  Distribution: 

HOB  (Meters): 

Fraction  of  Population: 

Precision  Error  Standard  Deviation,  Range 

Precision  Error  Standard  Deviation, 
Deflection 

MPI  Error  Standard  Deviation,  Range 
MPI  Error  Standard  Deviation,  Deflection 


850.  Ft/Sec 
30.° 

.95 

.3  5.  12. 

.03  .80  .17 

36.  meters 

6.  meters 
45.  meters 
27.  meters 


Surface  Height  of  Top  Edge  o 
Plate  or  Thickness  Inclination  to  Area  in  Surface  above  Ground 
Plate  Pair  Material  in  Inches  Horizontal  Inches?  in  Inches 
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.25  inches  hard  homogeneous  steel 
.20  inches  hard  homogeneous  steel 
.12  inches  7039  aluminum 


FRAGMENTATION  OATA  640.  G.032 
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FIGURE  03.  INPUT  DATA/  FRAGMENTATION 


TABLE  C-3.  RESULTS  OF  PROGRAM  SECTION  IV 
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PLATE 

PERFORATED 

P 

*I/H 

A 

yes 

.83 

B 

yes 

.62 
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no 
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no 
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no 
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no 
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yes 

.65 
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no 

_ 
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FIGURE  C-5.  OUTPUT  OF  PROGRAM  SECTION  IV  FOR  ONE  CF  THE  ARMOR  PLATES 


and  on  plate  H will  be  accounted  for  on  plates  B and  F and  on  plate  A, 
respectively.  Thus,  for  HOB  = .3  meters,  (HOB  assumed  to  be  half  the 
projectile  length  for  impact  functioning  rounds)  we  consider  plates  A, 

B,  F,  and  G,  and  for  HOB  = 5.  and  12.  meters,  we  consider  plates  A,  B, 
and  G.  It  should  be  obvious  that  plate  G will  sustain  little  damage. 

For  the  computations  of  program  sections  I and  II,  we  use  a 
four  inch  cell  at  HOB  = .3  meters  and  an  eight  inch  cell  otherwise.  We 
assign  a value  of  four  feet  to  the  blast  effects  sphere  radius.  We 
define  plate  A to  be  the  reference  plate  to  which  the  other  plates  will 
be  assembled;  therefore,  the  offsets  of  the  various  plates'  central 
cells  are  defined  with  respect  to  the  central  cell  of  plate  A (central 
cells  and  their  offsets  are  discussed  in  paragraphs  4.1  and  5.1).  Data 
for  the  individual  plates  of  the  structure  serving  as  the  model  of  our 
APC  are  at  Table  C-4.  The  input  data  deck  for  this  second  stage  of  the 
computations  is  at  Figure  C-6. 

Samples  of  the  output  from  program  section  I are  at  Figures 
C-7  through  C-9.  Table  C-5  contains  the  important  parameters  computed 
by  program  section  I. 

We  note  that  our  intuition  concerning  plate  G was  correct;  it 
has  sustained  no  damage.  We  note  further  that  for  HOB's  greater  than 
473.  inches  (55  eight  inch  cells  plus  33  inches  height  of  the  top  edge 
of  plate  G above  the  ground),  plate  G would  present  zero  surface  area  to 
the  burst  and  so  could  not  sustain  any  perforations.  Considerations  such 
as  these  can  be  important  when  machine  time  becomes  large.  The  source 
of  the  uncertainty  in  the  values  for  plate  B at  HOB  =12.  meters  is 
clear  from  examination  of  Figure  C-9;  perforations  might  be  sustained 
beyond  the  left  deflection  limit  of  the  matrix  used.  Any  error,  however, 
is  quite  small.  Finally,  all  these  values  refer  to  a half-infinite 
ground  plane. 

The  output  of  program  section  II  is  at  Figures  C-10  through 
C-12.  We  note  that,  due  to  the  possible  error  in  the  values  for  plate  B 
at  HOB  = 12.  meters,  the  12  meter  values  for  the  complete  structure  may 
incorporate  a small  error.  Again  these  values  refer  to  a half- infinite 
ground  plane. 

Figures  C-13  and  C-14  present  an  index  of  the  expected 
perforations  file  and  a tabulation  of  the  damage  function  file,  respec- 
tively. Although  not  produced  by  program  section  II,  these  are  presented 
here  because  the  input  data  card  which  caused  them  to  be  produced  was  a 
part  of  the  data  deck  shown  in  Figure  C-6.  The  first  several  logical 
records  (which  are  not  shown)  of  the  damage  function  file  pertain  to  the 
individual  plates  with  non-zero  probability  of  damage  (Table  C-5).  Damage 
function  file  logical  records  8,  9,  and  10  refer  to  the  complete  structure 
at  HOB  = .3,  5.,  and  12.  meters,  respectively.  Had  the  data  deck  of 
Figure  C-6  been  split  up  into  two  parts,  one  for  program  section  I and  the 
other  for  program  section  II,  then  an  examination  of  the  section  I 

Text  continues  on  page  193. 
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TABLE  C-5 . TABULATED  RESULTS  FROM  PROGRAM  SECTION  I 


HOB  in 

Lethal  Arf 

;a  Damage  Area 

in  Meters^ 

Average  Probability 

Plate 

Meters 

in  Meters* 

of  Damage 

A 

.3 

13.626 

13.626 

1 . 

5. 

9.323 

18.127 

.514 

12. 

8.930 

62.431 

.143 

B 

.3 

20.232 

20.232 

1. 

5. 

21.535 

95.835 

.225 

12. 

10.910+ 

171.521+ 

.064+ 

F 

.3 

11.613 

11.613 

1. 

12i} 

Not  computed; 

, known  to  be  0.  from 

section  IV  output 

G 

.3 

0. 

0. 

0. 

5. 

0. 

0. 

0. 

12. 

0. 

0. 

0. 
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THESE  DATA  ARE  ON  THE  DAMAGE  FUNCTION  FILE  (OFCN3A8  ) AS  RECORO  NUMBER  10 
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results  would  have  shown  that  plate  G could  have  been  eliminated  from 
the  section  II  deck,  since  it  contributed  nothing.  Although  very  little 
would  have  been  gained  in  this  example,  in  some  cases  this  sort  of  "one- 
step-at-a-time"  analysis  can  save  significant  amounts  of  machine  time. 

We  have  already  obtained  some  results  of  value  for  our  APC 
model:  lethal  areas,  average  damage  probabilities,  probability  of 
vulnerability  of  a single  APC.  Now  we  will  use  program  section  III  to 
obtain  measures  of  the  damage  sustained  by  the  APC  force  fired  upon  by 
the  howitzer  battery  (Figure  C-l).  Figures  C-10  through  C-12  provide 
contour  plots  for  the  damage  functions,  and  we  see  that  in  each  case, 
the  damage  function  is  "fully  nested"  (see  paragraph  6.1).  The  damage 
functions  input  to  program  section  III  are  constructed  from  those  output 
by  program  section  II  using  judgement  and  experience. 

In  this  example,  we  make  two  different  approximations  to  the 
damage  functions  output  by  program  section  II  in  order  to  give  some 
rough  idea  of  how  results  for  program  section  III  might  depend  upon  the 
input  damage  function.  The  input  data  deck  is  at  Figure  C- 1 5 . The 
first  approximation  results  in  "cookie  cutter"  damage  functions;  that 
is,  all  the  different  probability  of  damage  regions  are  combined  to 
obtain  one  region  with  a constant  probability  of  damage.  For  the  second 
approximation,  we  rectangul arize  (maintaining  normal izeticn)  the  various 
probability  of  damage  regions,  without  combining  any  of  them. 

Results  of  program  section  III  for  these  approximations  to 
the  damage  functions  are  at  Figures  C-l 6 and  C-l 7 . One  should  not 
assume  that  the  differences  in  results  for  different  approximations  to 
the  damage  functions  will  always  be  as  small  as  they  are  in  this 
example. 
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APPENDIX  D 


CONTROL  DECK  FOR  THE  EXAMPLE  OF  APPENDIX  C 


This  appendix  presents  and  explains  the  control  deck  used  to 
execute  the  program  in  the  example  of  Appendix  C.  Detailed  explanations 
of  the  commands  and  their  parameters  may  be  found  in  CDC  system  literature. 

Figure  D-l  presents  the  control  deck  used  to  obtain  the  first  set 
of  results,  the  program  section  IV  output.  Card  3 attaches  the  permanent 
file  on  which  the  program  is  maintained  in  relocatable  form.  Cards  4 and 
5 make  accessible  the  permanent  file  on  which  the  data  deck  of  Figures  C-3 
and  C-4  is  maintained.  These  data  could  be  read  in  on  cards,  but  they  are 
maintained  as  a file  for  convenience.  Card  7 causes  the  program  to  be 
loaded  and  executed  with  file  S4  replacing  the  system  input  file. 

Figure  D-2  presents  the  control  deck  used  to  obtain  the  second 
set  of  results,  the  output  of  program  sections  I and  II. 

Cards  4 and  5 make  accessible  the  permanent  file  on  which  the 
data  deck  of  Figures  C-3  and  C-6  is  maintained.  Again,  cards  could  be  used. 
Card  9 causes  the  program  to  be  loaded  and  executed  with  files  S12,  PMAT348, 
and  DFCN348  replacing  the  system  input,  perforation  matrix,  and  damage 
function  files,  respectively.  Card  10  causes  the  perforation  matrix  file 
to  be  catalogued  (made  permanent)  as  file  PMAT348  with  a password  whose 
future  absence  will  prevent  all  access  to  the  file  except  for  reading  it. 
Card  11  performs  the  same  function  for  the  damage  function  file,  under  the 
name  DFCN348.  If  an  additional  run  of  program  sections  I and/or  II  were 
desired,  with  additional  perforation  matrices  and  damage  functions  being 
written  on  these  same  files,  appropriate  ATTACH  or  GETPF  commands  would  be 
required  between  cards  7 and  9. 

Figure  D-3  presents  the  control  deck  used  to  obtain  the  third 
set  of  results,  the  output  of  program  section  III. 


ADDENDUM 


Paragraph  4.1.1,  page  8 , addresses  the  use  of  different 
linear  combinations  of  the  perforations  sustained  assuming: 

a.  All  fragments  within  a spray  leave  the  source  point  with 
that  spray's  median  velocity. 

b.  All  fragments  within  a spray  leave  the  source  point  with 
that  spray's  maximum  velocity. 

Subsequent  to  the  completion  of  this  report,  the  program  code 
was  modified  to  allow  the  linear  combination  to  be  varied  routinely,  as 
input  data.  The  coding  which  effects  this  is  shown  in  Figure  ADD-1. 

The  sequence  numbers  (columns  82-87)  indicate  placement  of  the  cards 
within  the  program  shown  in  Appendix  A. 

The  variable  linear  combination  coefficients  are  F(l)  and  F(2). 
F(l)  pertains  to  the  perforations  obtained  with  median  velocities  and 
F(2),  to  those  obtained  with  maximum  velocities.  F(l)  and  F(2)  are  read 
in  as  real  numbers  with  explicit  decimal  points  on  card  number  1 of  the 
BASIC  DATA  subdeck  (see  Figure  4-9),  F(l)  in  columns  71  thru  75  and 
F(2)  in  columns  76  thru  80. 
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